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Abstract Observations are reported of H2 IR emission in the S(l) v=l-0 line at 2.121 /im in the Orion Molecular 
Cloud, OMCl, using the GriF instrument on the Canada-France-Hawaii Telescope. GriF uses a combination of 
adaptive optics and Fabry-Perot interferometry, yielding a spatial resolution of 0.15" to 0.18" and a 
velocity discrimination as high as 1 kms~^. Thanks to the high spatial and velocity resolution of the GriF data, 193 
bright H2 emission regions can be identified in OMCl. The general characteristics of these features are described 
in terms of radial velocities, brightness and spatial displacement of maxima of velocity and brightness, the latter 
to yield the orientation of flows in the plane of the sky. Strong spatial correlation between velocity and bright H2 
emission is found and serves to identify many features as shocks. Important results are: (i) velocities of the excited 
gas illustrate the presence of a zone to the south of BN-IRc2 and Peak 1, and the west of Peak 2, where there is 
a powerful blue-shifted outflow with an average velocity of —18 kms~^. This is shown to be the NIR counterpart 
of an outflow identified in the radio from source I, a very young O-star. (ii) There is a band of weak velocity 
features (<5 kms~^) in Peak 1 which may share a common origin through an explosive event, in the BN-IRc2 
region, with the fast-moving fingers (or bullets) to the NW of OMCl. (in) A proportion of the flows are likely 
to represent sites of low mass star formation and several regions show multiple outflows, probably indicative of 
multiple star formation within OMCl. The high spatial and velocity resolution of the GriF data show these and 
other features in more detail than has previously been possible. 

Key words. ISM: individual objects: OMCl, source I - ISM: kinematics and dynamics - ISM: molecules - shocks 
waves 



1. Introduction 

The Orion Molecular Cloud (OMCl) (D =460 pc, 
iBallv et al.l 12000) is the closest highly active massive 
star-forming region. It has accordingly been studied over 
a large ran ge of wa velengths from the radio to the 



X-ray (e.g. JCenzel fc Stutz ki 1989: Wright et al.lll996t 
'Delll l200lt_JChrysostomou et all 120001 JCarmire et alJ 



o 



200rf 'Ferland"2001^. 'Feigelson et alJl200a iDoi et al.ll2002t 

O'DeU & Doi 2003: Tan 2004) and has become the type- 
site for studies of star formation. Many features are ac- 
cordingly very well-described, for example the detailed 
morphology of the HII region associated with 0^Ori C, the 
dominant sou rce of ionizing radiation in the Trapezium 
llO'Deli200ll) . Perhaps the principal significance of OMCl 
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and its vicinity is that it has been shown to harbour up- 
ward of 1000 young low m ass stars in the p resence of a 
number of high mass stars jHillenbrandlll99it . 

The present work, covering the area shown in Figs.^ 
and 121 concentrates upon infrared observations of molecu- 
lar hydrogen (H2 ) , and in particular upon data displaying 
the dynamics of OMCl. The gas in OMCl is highly struc- 
tured with many bright knots and there is ample evidence 
of much violent motion wit hin OMCl over a range of scales 
from 70 AU to 3 x 10"^ AU llOustafsson et al.l200fil). Earher 



work concerning gas dynamics (|£ !hrvsost omou et al.ll^ 
ISalas et alJIlQQa lOustafsson et al...2003.) showed the pres- 
ence of numerous gas flows with flow velocities varying 
from a few k ms~^ in relatively dense gas to severa l hun- 
dred kms-i l|Lee fc BurtonI 1200(1 IPoi et al.ll2002() . The 
latter are associated with perhaps the most prominent fea- 
ture in OMCl, the H2 "fing ers" (a l so kno wn as "bullets") 
first observed by lAUen fc BurtonI l)l993(l . The fingers lie 
principally to the NW of OMCl and they are the wake 
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of fast-movin g, diffuse clumps of gas moving radially out 
from OMCl llBurtonlll997tlLee fc Burtonll2n0f]l:lDoi et all 

l2002|l . We do not find any high velocity fingers in the ob- 
served field. This may be due to the limitations of the 
observations which confine us to velocities below about 
150 kms"i but may also be due to the technique of data 
reduction, as further described in sec. 12 

Two main categories of models have been proposed 
to explain the overall structure of OMCl. One is that 
structures in OMCl are caused by shock instabilities 
which form t he bright knots and bullets a s inlStone et alJ 
lll99,5 h and iMcCangh rean fc Mac Lowl l)l997l) . In this 
model Rayleigh- Taylor instabilities form in the wind from 
the central massive star (or stars) in OMCl when a faster 
wind overtakes the dense shell of material swept up by an 
earlier, slower wind. This results in a fragmentation of the 
shell into smaller knots some of which appears as bullets, 
forming fingers of emission, when moving outwards from 
the shell at high speed. The other gr oup of models (e.g. 
lAllen fc Burtonlll993t iDoi et alJlioO^ explain the bright 
knots as internal shocks in the cloud. These can have many 
origins: turbulent bipolar outflows from low-mass star for- 
mation or some form of outflow f rom the central mas- 
sive stars (BN, source I, source n llMenten fc R.eidlll99,'Tl 
iGreenhill et"al]l2004at IShuping et al.ll2004 ). Thefingers 
and bullets themselves in either model are typically inter- 
preted as the result of some major explosive event involv- 
ing one or more of the massive stars. The fingers may then 
have shot through Peak 1 rather than being formed there. 
The lack of fingers and bullets to the SE of OMCl could 
suggest a highly directional origin, but may equally well 
be explained by the emitting H2 being too deeply buried 
to be seen, or it has been evaporated by the action of the 
Trapezium stars; if you placed the bullets to the NW in 
the corresponding position to the SE, they would have 
come very close to the Trapezium. For further discussion 
of the fingers and bullets we refer the reader to f o r exam- 
plelAxon fc Tavloillll9 84'l: All en fc Burtonl lll99.'f):'BurtoE 




O'Dell et al.r(n997.) : .T.ee fc Burtonl VoOOil : lO'DeT 



Doi et all (j^ 



In this work we would favour the second group of mod- 
els mentioned above, in which a significant proportion of 
the excitation is internal to Peak 1 (see fig. [21 for nomen- 
clature). We also note there are no prominent fingers in 
Peak 2 and therefore Peak 2 must presumably have some 
internal excitation source. In separate work we have also 
studied the global natu re of the gas motions in OMCl 
i Gustafsson et alJl2006|) . In that work a statistical analy- 
sis is presented of the turbulent velocity structure in the 
Orion Molecular Cloud at scales ranging from 70 AU to 
3 X lO'^ AU, using the data of the present paper. The 
variance for the velocity field shows only an approximate 
power law dependence on lag, with deviations below 2000 
AU which are attributed to outflows and possibly disk 
structures, below 300 AU, associated with low mass star 
forma tion within OMCl . The major point is that the anal- 
ysis in lGustafsson et all l|20061 demonstrates the presence 
of one or more preferred scales in the excited gas. These 



preferred scales are most naturally attributed to star for- 
mation where deviations from power law behaviour is as- 
sociated with injection of energy at and below the star- 
forming scale. 

Energy injection into OMCl is evidently triggered 
within the molecular cloud at differ ent scales. Massiv e 
stars such as the BN object, a B-star 



Bzari et alJ Il? 

or radio source I (hereafter sirn p ly sou r ce I), a deeply 
buried 0-star ("Menten fc ReidI Il995t iCreenhill et all 

2004a; Shuping ct al. 2004), are associated with structure 
at a scale of ~0.1 pc (20,000 AU). Low mass star forma- 
tion gives rise to flows and associated structure at scales 
of hundreds to several thousand AU. 

The region studied in the present work is shown in Fig. 
1 and lies NW of the Trapezium stars. The target area cov- 
ers only the inner part of OMCl, and as we have noted 
largely excludes for example bullets and fingers. Thus the 
material in the present study encompasses the gas volume 
that tends to be subject to denser, slower outflows, associ- 
ated with both massive star formation and low mass star 
formation. 

High spatial resolution K-band (2-2.5 ^m) images 
of IR emission from in OMCl u si ng the HST 
llStolovy et a,].l 1799^ IChen et a,lJ ISchnltz et alJ 

Il999t iDoi et alJ I2OO2I). the Canada-Fr ance-Hawaii 
Telesc ope (CFHT) (IGustafsson et al. 2003^ the ESO 
3.6m Vannier et a'ill2'o01t iKristensen et alJ l2003(l . the 
VLT ( Lacombe et al.l I2OO4I) aiid other telescopes (e.g. 
IMcCaughrean fc Mac Lowlll997t ISchild et alJll997l) have 
provided the most recent evidence for the presence of fast 
flows and shocks at small scales in OMCl. 

The recent paper of IGustafsson et al.1 l)2003|) . here- 
after Paper I, provided an initial report of H2 emission 
data for OMCl obtaiiied us ing the so-called "GriF" in- 
strument l(Clenet et al.ll2002l) on the CFHT in December 
2000. This instrument combines Fabry-Perot interferome- 
try with the PUEO adaptive optics system on the CFHT 
i Rjgaut et al.lll998j) . These data, containing information 
on both morphology and velocity, showed striking evi- 
dence for the presence of shocks within OMCl. This was 
provided through the clear association of gas motions of 
tens of kms~^ with bright IR emission of H2. The nov- 
elty of these data is contained in the association of a 
very high spatial resolution of 0.15" (70 AU) to 0.18" 
with velocity discrimination as good as 1 kms^^ (3cr) 
in regions of high emission brightness, obtained with a 
Fabry-Perot of relatively low resolution of 150 kms~^. 
These figures may be compared wi th the high inherent 
spectr al resolution of 14 kms~^ in Chrvsostomoii et alJ 
l)l997l) and 24 kms'^ in ISalas et all l|l999l) . where data 
were recorded with a spatial resolution which lay be- 
tween 1.5" and 2.0". In the present work the GriF data 
are reported in detail. Whereas in Paper I 19 regions 
were identified for study, 193 regions are included here. 
This affords a much broader perspective on the nature 
of flows within OMCl and represents a detailed attempt 
at extracting the wealth of information contained within 
the GriF observations. New features are revealed through 
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Figure 1. The observed area relative to the Trapezium 
stars (large black dots). The grey scale image show emis- 
sion brightness in the H2 v=l-0 S(l)line. Axes are labelled 
in arcseconds. The star TCC0016, at coordinates 0" E, 
0" N and marked with a small black dot, is used as a 
reference for coordinates throughout this work. 




Figure 2. The regions analysed in this work outlined as 
small overlapping white boxes. The gray-scale shows emis- 
sion in the H2 v=l-0 S(l) emission line. Also shown as 
large wh ite boxes is the nami ng convention for Peak 1, 
Peak 2 ijBeckwith et al.lll978|) . and Region B, the lat- 
ter bounded by -2" --40" E, -18" -17" N, relative to 
TCC0016. Axes are labelled in arcseconds. 



this extended study of the GriF results. In particular, a 
massive and very energetic blue-shifted outflow is iden- 
tified in Region B between Peak 1 and Peak 2 (fig. 
IChrvsostomou et"al] l|l997|) . working in K-band at a spa- 
tial resolution of 1.5" found evidence of clumps of gas dom- 



inated by blue-shifted motions in the same zone. As de- 
scribed in sect. El there is a considerable independent body 
of evidence from radio data for an outflow in this zone 
fcenzel et al 1981; Genzel fc Stutzki 1989: Wright et alJ 
Il996t iGreenhill et al.ll2004al lOoelemanet aiJl2004() . The 
outflow detected here in excited H2 is characterised in 
some detail and it is shown that the outflow is the IR 
signature of the outflow identified through radio observa- 
tions. 

On a smaller scale, numerous regions of bright H2 emis- 
sion are identified in Peak 1, the NW region, which show 
only weak radial velocity components with no clear veloc- 
ity structure. In addition there are several regions on the 
1000 AU scale with multiple outflows. If, as we discuss 
in some detail here, there are protostars associated with 
some of these latter zones, then our data present a picture 
of star formation that is more complex than the simple 
disk-bipolar outflow model of single star formation. 

In sect. 13 the observations and data reduction, de- 
scribed in detail in Paper I, are briefly reviewed. Sect. |31 
illustrates and collates the data for all 193 regions in terms 
of radial velocities of flows, brightness of emission associ- 
ated with these flows and the displacement in the plane of 
the sky between the maximum in emission brightness and 
the maximum in local radial velocity. The mechanism of 
H2 excitation is also briefly discussed. In sect.0I the data 
are presented in order to give a global view of the spatial 
distribution of radial flow velocities and the characteris- 
tics of various zones. Among major points that emerge are 
(i) the presence of the outflow mentioned above between 
Peaks 1 and 2, in Region B (Fig. [5J, (ii) the accumula- 
tion of zones with high II2 brightness but with <5 kms~^ 
radial motion in a region passing diagonally across Peak 
1, (Hi) the greater tendency for shocked regions to group 
more closely in Peak 2 than in Peak 1 and (iv) the fact 
that maxima in brightness of emission and in velocity are 
in general slightly displaced. In sect. [3 a short qualitative 
account of the structure of shocks is followed by the use 
of data to establish the orientations of shocks within the 
plane of the sky. Then follows a consideration of possi- 
ble sites of current low mass star formation in OMCl. In 
sect. El the origin of the outflow in Region B and its 

spatial association with continuum sources in the IR 
and radio is discussed in detail. 



2. Observations and data reduction 

Observations were performed of the regions shown in 
Figs. [Hand El on December 5th 2000 at the 3.6m CFHT , 
using the then new instrument GriF l|Glenet et alJl20o3) . 
This instrument combines the high spatial resolution of 
the PUEO adaptive optics with spectral resolution pro- 
vided by a Queensgate ET50WF Fabry-Perot interferom- 
eter (FP), giving a data-cube showing brightness in each 
position (pixel) for a series of wavelengths. The total ob- 
served region in OMCl is a ~ 1' x 1' area centered approx- 
imately on the BN-object, as shown in Fig. [3 
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Observations were made in the K-band using the 
H2 v=l-0 S(l) e mission Une, which has a rest wavelength 
of 2.1212544 ^m (|Bragg et al.ll982l) . The line was scanned 
using the FP with a step size of ~4.5 x 10""* fxm to obtain 
images showing the brightness at a series of wavelengths 
on both sides of the rest wavelength. A H2 v=l-0 S(l) in- 
terference filter of central wavelength 2.122 /im and band- 
width of 0.02 ij,m was inserted between the FP and the de- 
tector to prevent the superposition of different FP orders. 
The spatial resolution of the data was 0.15", estimated 
from the PSF of stars in the field. The pixel scale was 
0.035". The inherent spectral resolution of the FP was 
150 kms^^. To reduce the level of noise, the data have 
been smoothed with a moving boxcar average over 3x3 
pixels, resulting in a 20% loss of spatial resolution, that 
is, the final resolution of the data is 0.18". For a given po- 
sition in the plane of the sky the data give the count rate 
for each o bserved wavelength . Thi s is the n fitted with a 
lorentzian l|Clenet et al.ll2002l: IChrv sosto mou et alJll997l 
and Paper I) to find the maximum count rate and the cor- 
responding wavelength. This gives the velocity of the gas 
at that position. Due to high signal-to-noise very precise 
fits can be made. In bright regions the uncertainty in the 
fitted velocity is ±1 kms~^ or less at 3a. In the dimmest 
regions examined in this work (with a very few exceptions) 
the uncertainty in the peak in the lorentzian, used to de- 
termine flow velocities (see sect. 13. it . is <2 kms""'^ while 
velocities of the gas surrounding regions of flow typically 
have uncertainties of 3-4 kms~^. A detailed account of 
the error s associated with velocity determination may be 
found in iGustafsson et al . ( 2006) where the uncertainty 
associated with velocity determination and the counts per 
pixel per exposure are explicitly related through an em- 
pirical formula. 

Where weaker high velocity emission is seen spatially 
coincident with stronger low velocity emission, we record 
only the velocity of the brighter slower gas. This is however 
a rare occurence and we are aware of only two cases where 
this has been encountered: (i) th e emission feature ar ound 
position -24.2" E 6.5" N where IStolovv et alJ ^^1998^ find 
a weak high velocity component together with the stronger 
low velocity feature which we find and (ii) at —19.7" E 
16.2" N where we observe the object HII208 at a radial 
velocity of —19 kms"'^. HH208 has however been identi- 
fied as a bullet, on the basis of associated Fell emission, 
with a radial velocity between —120 and —180 kms~^, 
lAxon fc Tavloil lll984l) : lO'Dell et all ^ll997^ . 

Velocity calibration was carried out by assigning to 
the mean velocity of all data a value of 12±6 kms~^ in 
the local standard of rest ( y;sr), consistent w ith data in 
IChrvsostomou et alJ l)l997|) : ISalas et al.ll)l999l) and lO'Delll 

Velocity calibration is however not material to this 
work, as we are concerned with relative velocities over 
small distances of typically a few arcseconds. For further 
details on observations and data reduction see Paper I and 
references therein. 



Throughout this paper all positions are given in arcsec- 
onds east and north relative to the star TCCG016, which 
itself is situated at 05''35"n4!91, -05°22'39"31 (J2000). 

3. Data analysis 

3.1. Identification and characterization of emitting 
zones 

Data pertaining to the 193 regions analyzed here, are col- 
lated in tables rO - rOl All flows are designated by num- 
ber in these tables. The tables each refer to one of the 
three regions into which the full region has been subdi- 
vided, namely Peak 1, Peak 2 and Region B. 

Features were identified on the basis of their brightness 
and our concomitant ability to locate region boundaries. 
Thus all regions show a well-defined localized increase 
in brightness compared to the surrounding gas. Each of 
these regions has been examined for evidence of a local 
gas flow. This was done by making 3D-images of each re- 
gion, where these images show radial velocity and emis- 
sion brightness as a function of the position in the plane 
of the sky. Two illustrative examples of regions exhibit- 
ing flows are shown in Figs. and b, where the x- and 
y-axes give the position in the plane of the sky and the 
z-axis shows the radial velocity of the gas. A colour code 
associated with every image shows the observed brightness 
as the number of counts per exposure of 400 s. This can 
be converted to brightness through the us e of the value 
of 3.0 ±0.15 X IQ-^ Wm-^sr-i found in IVannier et alJ 
l|200l|) for the brightest region observed here.^ Using this 
value it is found that a count of 100 over the 400 s exposure 
corresponds to ^^3.5 x 10^^ Wm^-^sr"^. This conversion 
may be applied to all regions (and all figures in this work) 
with the proviso that systematic errors may arise through 
differing conditions of airmass and of the atmosphere pre- 
vailing at the times of observation. 

Fig. clearly shows a very localized increase in the 
radial velocity^ and the brightness. The typical extent of 
the areas which show this phenomenon of strong localised 
brightness is 0.5" -1" in all the 193 regions identified in 
the data. All regions are therefore well resolved, given a 
spatial resolution of 0.18". 

Data for each region have been characterised in terms 
of the following properties. 

(i) the maximum velocity relative to the surrounding gas, 

(ii) the maximum brightness, 

(Hi) the displacement in the plane of the sky between the 
maxima in velocity and brightness. 

(i) The maximum velocity of the flow relative to 
the surrounding gas is found by taking the maximum 
velocity recorded in the region of interest and subtracting 
the average velocity of the surrounding gas. The latter 

^ The brightest region is in Peak 2 at coordinates 15.82" E, 
-1.61" N, see tableEH flow no. 2-10. 

In the remainder of this article, velocity should be taken to 
mean the radial velocity, unless otherwise specified. 
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Figures, a Flow no. 1-61: see table lA.ll The vertical axis shows v/sr in kms^^. Colours denote brightness as shown 
on the colour bar above, in counts per 400 sec. The plane of the sky is shown in arcseconds, relative to TCC0016. b 
Flows no. B-44, B-45 and B-46: see table IXTsl Otherwise as in a. c Flow no. 1-74: see table IXtI This is an example 
of a flow which shows no clear velocity peak exceeding 5 kms^^. 



may be estimated with a typical uncertainty of 3-4 kms~^ 
(see sect. The resulting flow velocity contains this 
uncertainty. Thus these velocities are less accurately 
determined than relative velocities between adjacent 
bright regions. Using Fig. as an example, the velocity 
of the surrounding gas is 8±3 kms^^. The tip of the 
velocity peak is at — 24±1 kms~^ and thus the flow 
velocity relatively to the surroundings is — 32±4 kms~^ 
(3a). 

In 64 out of the 193 regions there is an increase in 
brightness similar to that found, for example, in Figs. Eli 
and b, but with no corresponding velocity structure. 
An example is shown in Fig. These features will be 



referred to as features with weak radial velocity compo- 
nents (WRVs). Fig. int, which is a typical such object, 
demonstrates that velocity variations are of the order 
of less than 5 kms~^, with no clear velocity structure. 
This is therefore used as a definition of WRVs, which 
are defined as those objects with a localized maximum 
in brightness but which show no velocity excursions 
from the surrounding gas greater than 5 kms~^. WRVs 
according to their brightness and extent may arise from 
flows close to the plane of the sky. However, some of these 
features could also be bright photodissociation regions 
(see sect. I3.2|l . In this connection data for the great 
majority of WRVs are of sufhciently high signal-to-noise 
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that variations in velocity are real and do not simply 
arise from uncertainties in velocity determination. 

(ii) The maximum brightness is the highest recorded 
brightness within the confines of any region after the 
3x3 boxcar smoothing mentioned in sect. El 

(in) Spatial displacements of brightness and veloc- 
ity features are listed in Tables lO - EH There exists 
a strong spatial correlation between the local peaks in 
brightness and velocity: the distance in the plane of the 
sky between the two peaks is typically less than the 
extent of either peak. That is, they overly each other in 
the great majority of cases, with the maximum in velocity 
lying within the area where the emission caused by the 
flow is more than than twice the background emission. 
Figs. ISt and b illustrate this general property that 
features of velocity and brightness are closely associated 
in the plane of the sky. The use of these displacement 
data to determine the orientation of gas motion in the 
plane of the sky is discussed in sect. 15.11 

3.2. The Origin of the H2 emission 

The mechanism by which the H2 emission is formed 
is brie fly discussed here. This has already been cq nsid- 
ered in lVannier et alJ ll200lh.lKristensen et all l|2n0,'^ and 
Paper I. 'Vannier et al. (2001) described in detail, using 
the shock models of Wilgenbus et al. ( 2000), how the 
brightest regions in the Peak 2 could arise from mag- 
netic (C-type) shock compression of dense gas at pre- 
shock densities around 10^ cm~'^ by shocks travelling at 
30 kms^^. Using sh ock models and limit ing their inves- 
tigations to Peak 2, IVannier et alJ l)20 0l') sought to find 
if quantitative evidence could be established that shocks 
have collected enough gas and compressed it sufficiently 
to overcome the Jeans criteria. This work was therefore 
directed towards discovering if shock-induced star for- 
mation is an active mechanism. Later VLT observations 
have in fact shown the presence of a 140 AU disk in L- 
band around th e position of b r ightes t emission (15.8" E 
-1.6" N) where fVannier et alJ l)200lf) concentrated their 
attention. Therefore a protostar would appear already to 
have formed at this position and the emission observed 
was most likely due to outfiow from a YSO rather than 
the much earlier stage of shock induced compression en- 
visaged in IVannicr et al. (2001). While other candidates 
in Peak 2 may be less compelling in terms of emission 
brightness and the amount of material within a restricted 
volume, the mechanism of shock-induced star formation 
could nevertheless be active elsewhere Peak 1 and Peak 2. 

The shock-induced mechanism leads of course to star 
formation at some time 10^~^ years in the future, assum- 
ing collapse on a free-fall timescale. In this connection, 
according to the latest shock models Ip^^ Bourlot et al.' 
l2002t iFlower et all l2003t iKristensen etalj .2003. .2006^) . 
the brightest emission (>10~^ Wm~^sr~^, >300 counts 



per 400 seconds) may readily be generated by magnetic 
shocks of typically 25 k ms~^ impinging on d ense gas 
(>10^ cm~^), much as in IVannier et all (|200ll) . Shocks 
leading to the forr nation of the fing ers occurred recently 
(~ lO'^ years ago, iDoi et all (|2002fl ') and the associated 
shock in dense gas in inner regions of OMCl, such as Peaks 
1 and 2 may lead to further star formation, again on the 
iC^^ year timescale. 

Both IKristensen et al ] (Eooi and Paper I addressed 
the question of the relative importance of shocks and 
photodissociation regions (PDRs) in yielding the ob- 
served H2 emission brightness. 6'^Ori C generates a 
UV field, at around 100 nm, more than 10^ times 
greater than the standard interstellar field. Detailed mod- 
els for H2 emission have been devised, for example in 
IStorzer fc HoUenbach (199B,), for UV fields of this inten- 
sity falling u pon dense gas, including advection and photo- 
evaporation teennev & 0'Del]lll99a) . These models show 
that the maximum contribution to emission brightness in 
the H2 v=l-0 S(l) line from PDRs docs not exceed a few 
times 10^^ Wm~^sr~^, given a line-of-sight normal to 
the PDR. For this geometry, the PDR contribution to 
the brightness is therefore no more than 10% to 15% of 
the maximum observed here. The conclusion is that the 
brightness due to photoexcitation is a generally minor, 
though non-negligible contributor to the H2 emission ob- 
served here. This issue is discussed in more detail, w ith 
reference to specific regions, in Kristcns en et al.l ((20031) . 

Possible caveats to the above could be provided by 
young highly obscured 0-stars around which the UV field 
would be even higher than the standard Orion value and 
the H2 would be of very high density. However the re- 
sulting H2 emission would also be highly obscured from 
our view. Low-mass YSOs of class 2 or 3 may also con- 
tribute with 100 nm emission, but a comparis o n wit h the 
catalogues of stellar sources in iMuen ch et al.l (|2002ll and 
iLada et al. (20Q4) (see sect. 15.21 and fig. (TTjl shows that 
stellar objects of class or I are only found near strong 
H2V=l-0 S(l)emission in a few cases. Note that with a 
few exceptions all of the areas examined in this work have 
a maximum brightness which is 25% or more of the ob- 
served maximum brightness, and all regions showing clear 
evidence of flows show brightness above the 15% limit 
mentioned above for PDR contributions. Thus PDR mech- 
anisms are unable to reproduce the observed brightness 
of H2 emission. This issue is re-examined briefly in the 
discussion of weak radial velocity components (WRVs) in 
sect. 15.31 where other geometries are considered in which 
PDR activity may play a more significant role. 

4. Results 

To illustrate the nature of the data collated in Tables fOI - 
IA.3I all features have been plotted in Fig. ^ where each 
circle marks an individual feature. The radius of the circle 
is proportional to the maximum brightness, and the colour 
shows whether features are red-shifted (red circles), blue- 




Figure 4. H2 v=l-0 S(l) emission at 2.121 fim (gray scale). Circles show all 193 features identified in the present work in this line. The centre of each circle 
marks the position of maximum brightness of a feature, and the radius of the circle is proportional to the maximum brightness. The colour of a circle shows 
whether a flow is red-shifted (red circles), blue-shifted (blue circles) or with no clear velocity feature, WRVs - see text, (yellow circles). The position angle of 
the flow, if determined, see sect. 15. l1 is shown by a black line in the circle. The star diagrams in the upper left hand corner show all position angles in a given 
region. The length of the lines in each star diagram corresponds to the displacement between maximum brightness and maximum velocity. Points P and Q on 
the main figure axes define the orientations of the line of WRVs across Peak 1 (see sect. 14.31 and sect. l5.3|l . 
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shifted (blue circles) or if WRVs were found to be present 
(yellow circles). 



4.1. Flow velocities 

Figure shows the number of flows for a given velocity 
for Peaks 1 and 2 and Region B. Several properties of the 
flows are apparent: 

(i) The fraction of red-shifted flows is similar in Peaks 1 
and 2, with 33% in Peak 1 and 45% in Peak 2. There are 
thus overall fewer red-shifted than blue-shifted flows. The 
detection of bright H2 emission is strongly affected by the 
local dust optical depth. In dense regions associated with 
star formation the ne ar-IR may be strongly obscured 
l)R,osentha,l et a,1. 1 1200(1) . Since the densities may exceed 
10'' cm~'^, one hundred to a few hundred AU depth 
of material is sufficient to obscure H2 emission. Thus 
material is preferentially observed emerging from the 
obscuring dust, rather than retreating into it, that is, blue 
shifted flows are expected to be detected more readily 
than red shifted, as found. 

(ii) Compared to Peak 2, Peak 1 has a large over- 
representation of WRVs, as evidenced by Fig. [S^ and as 
is clear from Fig.^J WRVs in Peak 1 may have a different 
origin than other features in Peaks 1 and 2, arising either 
from a concerted outflow associated with the BN-IRc2 
region and/or with a possible PDR contribution. This is 
discussed in sect. 14. 31 and sect. 15.31 

(Hi) the total number of flows in Peak I and Peak 
2 is 40 in each case, excluding WRVs since these may 
represent in general a different class of object. Since 
Peak 2 is less than half the size of Peak 1, the density of 
flows is correspondingly larger in Peak 2. 

(iv) In Region B all flows are blue-shifted and the 
velocities are distributed around an average value of 

— 18 kms~^ with a spread of 8 kms^^. 

(v) For both Peak 1 and Peak 2 the red- and blue- 
shifted flows appear clumped around -1-10 kms~^ and 
— 10 kms~^. Recall that velocities are relative to the visr- 
Thus this symmetric distribution is the result we would 
expect if the origin of the flows is internal to OMCl, 
be it turbulent motion or outflows from protostars. The 
velocity dispersion of p rotostar s is on l y a few kms~^ 
in th e Orion Nebula l|0'Del]l l200l[ Ivan Altena et all 
Il988() . Hence if H2 emission arises to some extent from 
outflows from protostars, one would expect to observe 
a distribution of velocities which reflects the range of 
outflow velocities, since these greatly exceed relative 
protostar velocities. Thus data for peaks 1 and 2 in fig. 
are consistent with the presence of protostellar outflows. 



(vi) Flow velocities lie between a few kms ^ and 



an upper limit of '^35-40 kms^^ for the radial velocity. 
Again, on the basis that some of the flows may be associ- 
ated with protostellar outflows, these figures agree with 
results from^avis et al. ( 2001) where fiow velocities of 5 
- 40 kms~^ are recorded for a sample of 9 protostellar 
zones. The up per limit is consistent with results in 
iLe Bourlot eTa l. (200^ for the critical velocity at which 
a C-type shock collapses to J-type, in gas of pre-shock 
density 10^-10^ cm^^. 

(vii) Fig. |S1 shows Hubble diagrams for each of the 
3 regions. Peak 1, Peak 2 and Region B, plotted as 
the radial velocities versus distance from source I. The 
dotted line in fig. represents an age of 1000 years 
co rrespondirig to t he age of the finger system as found 
bv lDoi et alJ l)200^ . There is no evidence of Hubble- type 
behaviour in any of the observed regions. Some of the 
data in fig. |^ may contain close-in representatives of the 
finger system, the great majority of whose members lie 
outside our field as noted earlier. It is clear however that 
the data provide no evidence that the observed shocks are 
part of the concerted outflow from a common explosive 
event with which the flngers are associated. One proviso is 
that Fig. El show only radial velocities. Thus the lack of a 
distance- velocity correlation does not in itself exclude an 
origin of the observed emission in some common explosive 
event. Thus we may observe random azimuthal velocities 
associated with flow in the plan of the sky, observed here 
as random radial velocities. In fact we show, by appealing 
to orientation angles of flows in sect. 15.1.21 that there is 
some directed flow in Peak 1. We should also recollect 
the present discussion excludes WRVs. In sect. 15.31 we 
find that there may be a common origin of the WRVs in 
Peak 1 and the bullets, again representative of a directed 
fiow to the NW. 



4.2. Displacement of velocity and H2 emission 
brightness 

Fig. shows the distribution of displacements between 
associated maxima in brightness and in velocity for each 
region. We find that for 2/3 of the flows the displacement 
between maximum velocity and maximum brightness is 
<0.6", and less than 10% show a displacement of >1.0". 
This should be compared to the typical extent of a peak 
in velocity or brightness which is of the order of 0.5"- 1". 
Data in Fig. [7| quantify the earlier statement that a 
strong correlation exists between the positions of maxima 
in brightness and velocity (sect.|OJ, implying that these 
are connected phenomena. 



4.3. Spatial distribution of flows with weak radial 
velocity components (WRVs) 

Several features are apparent in the spatial distribution 
of WRVs. 
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Figure 5. Histograms showing the number of flows with a given radial velocity relative to the velocity of the immediate 
surroundings in each of the three main regions. Velocities are binned over 4 kms~^, except for WRVs which are 
inherently binned over 10 kms^^(see sect. Weak radial velocity features (WRVs) are shown as gray-shaded bars. 



(i) WRVs comprise only a small proportion of the 
total number of flows in Peak 2 while in Peak 1 they 
outnumber all of the flows identified. 

(ii) If the assumption is made that the WRVs arise 
from shock excitation involving clumps of gas moving 
in or very close to the plane of the sky, then there is 
an inhomogeneous distribution of motion within Peak 1. 
This may be shown as follows. 

If the general flow pattern in Peak 1 were homogeneous 
in three dimensions, then an estimate may be made of the 
number of flows that are expected to fit into the cate- 
gory of WRVs. Given that typical shock speeds in OMCl 



are~30 kms~^ (as suggested by results in lVannier et alJ 
1200 ih . then the limit imposed on the radial velocity to 
define a flow as a WRV, ^radial < 5 kms~^, sect. 13.11 cor- 
responds to flows moving within approximately ±10° of 
the plane of the sky. This is illustrated in Fig. |2| For a 
homogeneous distribution of flow direction, the number of 
flows designated as WRVs relative to the total number of 
flows should equal the relative area of a band on a sphere 
around the equator from —10° to -1-10° (= 7r/18 rad.): 



no. of flows ill the 
plane of the sky 



J-n/18 cos 



t/18 



total no. of flows 4tt ^ ^ 

In Peak 2 a total of 50 flows (including WRVs) are found 
and thus in a model of homogeneously directed flows, 
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Distance from Source I (orcseconds) 

Figure 6. Hubble diagrams showing velocity versus distance from Source I for each of the 3 regions. In a, the dotted 
line co rresponds to an outflow of age 1000 years corresponding to the age of the OMCl fingers according to lPoi et alJ 
ll2002l) . Note that the WRVs are not included in these diagrams since these have no measureable radial velocities. 
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Figure 8. Illustrating the ratio of the number of flows 
showing less than 5 kms~^ radial velocity (WRVs) to the 
total number of flows: see sect. 14.31 and Eq. 1. 

8-9 WRVs are expected. This agrees well with the 10 
examples observed. The same model applied to Peak 1 
suggests that ~15 WRVs should be seen. However 47 
WRVs are observed in Peak 1. The above results do 
not change significantly if a typical shock speed of 25 
or 20 kms^^ is assumed instead of 30 kms^^. Thus the 
flows do not form a spatially homogeneous sample in 
Peak 1 but lie preferentially in the plane of the sky and 
appear disproportionately as WRVs. 

(Hi) The data in Fig. 0] show not only that there is 
a considerable excess of WRVs in Peak 1 but also that 
there is the appearance of a concentration of WRVs 
along a direction on the eastern side of Peak 1 from 



-8" E, 47" N to -31" E, 21" N along the hue connecting 
points P and Q in Fig. |3| This hue (or sht) is --30" long 
but only a few arcseconds wide. To assess whether this 
is a genuine feature of the data or merely a trick of the 
eye, a slit of length 30" and width 2.5" has been placed 
at all positions and for each position at all angles over 
the data for Peak 1 in Fig. 21 In each case the number of 
WRVs has been counted within the slit. Figure |51 shows 
the number of WRVs in the slit at the optimum angle in 
the sky to pick up the maximum number of WRVs, that 
is, 140° measured clockwise from N, as a function of the 
perpendicular distance of the slit from TCC0016. This 
corresponds to the slit being placed along lines parallel 
to the line connecting points P and Q in Fig. 0] The 
presence of a very clear maximum maximorum indicates 
that there is a genuine concentration of WRVs at the 
position designated. There is thus some kind of "front" of 
emission features at this position. 

(iv) A further property distinguishes the WRVs from 
other features in the field within this concentration or 
front just identified. An investigation of the morphology 
of the H2 emission in these features reveals that '---2/3 of 
these objects show a very much more broken up morphol- 
ogy than that typically seen in features elsewhere in Peak 
1 and Peak 2, for which emission is clearly associated with 
radial motion. Moreover, even in other WRVs, outside of 
this concentration, such a broken up morphology is only 
seen in 10-15% of the cases. The broken-up morphology 
associated with the concentration of WRVs is illustrated 
in Fig. llUb . This may be compared to the relatively 
uniform increase towards a single maximum seen in most 
other WRV objects such as the WRV shown in Fig. llUb . 
The origin of brokcn-up structure is briefly considered in 
sect. 15.31 
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Figure 7. Numbers of flows with a given displacement between the maxima in brightness and velocity for each of the 
three main regions. 



4.4. Grouping of flows 

In the present section all flows in Fig.^are considered save 
WRVs. Fig.^]shows how many blue- and red-shifted flows 
lie within a distance of 1200 AU (2.6") from any given 
point, where each flow is represented by a circle of radius 
1200 AU. A radius of 1200 AU was chosen because a re- 
gion of 2000-2500 AU is a typical scale for a s tar forming 
core in the early stages of star forma tion, e.g. iBate et all 
l)2003|) : iDelgado-Donate et alJ t004\ and this is also the 
scale where lGustafesoirer^Lri|2006|) find deviations in the 
turbulent velocity structure, showing the presence of a pre- 
ferred scale. Fig. ^2 illustrates whether the flows have a 
tendency to be grouped together or are spaced equally 
across the entire area. Peak 1 and Peak 2 show different 
character. In Peak 1 the flows are distributed rather evenly 



across the entire area, with only one position showing a 
large number of flows in a small area, around the position 
—22" E, 36" N. In Peak 2, flows are more grouped with 
some areas showing 6-8 flows within an area of ^2500 AU 
in diameter around positions 1" E, 5" N and 9" E, 1" N. 
Taken together with the observation, reported in sect. 14. II 
that the spatial density of flows is approximately twice as 
high in Peak 2 as in Peak 1, this points towards Peak 2 
being a more active region compared to Peak 1: see also 
sect.lO 

4.5. Further characteristics of the outflow in Region B 

As noted, all flows are blue-shifted in Region B and the 
mean of the velocities is significantly greater than in 
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Figure 9. The number of WRVs in a slit placed across 
Peak 1 at an angle of 140° measured clockwise from N as 
a function of perpendicular distance between the slit and 
TCC0016. The position of the peak at 37" corresponds to 
a slit placed along the line connecting points P and Q in 
Fig.H 

Peaks 1 and 2. Moreover Region B contains, among the 49 
objects in total, five clear examples of bow shocks, based 
on the morphology of the H2 emission. An example is 
displayed in Fig. 1121 the left-hand panel of which shows 
a NAOS- CONICA adaptive o ptics image obtained with 
the VLT ijLacombe et al.ll2nn4l) . lacking velocity data but 
showing the morphology in graphic detail, and the right 
hand side showing the GriF results, including velocity data 
as contours. 

The structure in Region B appears at first to con- 
sist of disconnected fragments of emission. However on 
closer examination it is found that there is widespread 
weak emission connecting some of the prominent bright 
features shown in Fig. ^ An example is shown in Fig. 1131 
A ridge of very weak emission, at the level of 30-60 counts 
per 400 seconds, corresponding to velocity uncertainties 
of 3 - 6 kms^^, can be identified connecting three bow- 
shocks showing blue shifted velocities between —10 and 
—30 kms^^. Thus the outflow region represents a con- 
certed bulk motion, though diffuse in places. 

5. Discussion 

Gas flows within OMCl are of two types. The first is 
large scale flows. Two such flows may be identified in 
OMCl. There w as a large-scale outburst ^1000 years ago 
()noi et a,iJl2nnl roughly in the plane of the sky, which 
has lead to the subsequent formation of fast moving fin- 
gers aiid__buJjBtsjn th*^ to the NW 
e.g. iLee fc Burtm] l|2000l) : IPoi et all l|2002l) . The WRVs 
concentrated in Peak 1 may be part of this scenario, as we 
discuss in sect. 15.31 The flow identified here in Region B 
constitutes the other large scale feature. By contrast, this 
is largely out of the plane of the sky. The relationship be- 
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Figure 10. a Flow no. 1-57: see Table IXm A smoothly 
structured WRV situated in Peak 1 but not in the flow 
front described in the text: sect. 14.31 b Flow no. 1-7: see 
Table IXTl A highly structured WRV feature lying in the 
flow front. Colours denote brightness as shown on the 
colour bar above, in counts per 400 sec. 

tween these two outflows is not presently clear. The out- 
flow in Region B is discussed in detail in sect.lHI 

The second type of feature is small-scale flows, asso- 
ciated with local low mass star formation. Flows in this 
case, taken as a body, would be expected to show no over- 
all preferential direction of motion within the plane of the 
sky, with each site or group of sites acting independently 
of any other. This point is examined in sect. 15. l1 below. 

5.1. The physical nature of shocks in OMCl and their 
orientation in the plane of the sky 

This section considers how data for the displacement of the 
maxima in velocity and brightness in any chosen feature, 
recorded in the final column of Tables lA.ll - IA.3I may 
be used for the determination of the orientation of shocks 
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Figure 11. Spatial density of emission features associated with flows, showing the number of flows within ~1200 AU 
of any given point. The colour scale illustrates the number of overlapping features at any point. See sect. 14.41 and 
sect. 15.21 Axes are labelled in arcseconds. 



in the plane of the sky. To enable this interpretation, it is 
necessary to consider the brightness and velocity structure 
of C-type shocks and the spatial relationship between the 
regions of maximum velocity and maximum brightness in 
H2 emission. 



5.1.1. A brief overview of shock structure 

The following gives a simplifled description of the general 
events at play in a shock. This description is only applica- 
ble to features in OMCl in broad outline but presents 
a physical picture which should prove hclpfulin inter- 
preting our observations. We envisage the shock form- 
ing as follows: in a time dependent picture, starting at 
time zero, before the shock has been created, a jet of 
gas, from hereon for simplicity considered continuous, is 
ejected from a protostar and collides at a super- Alfvenic 
relative velocity with a dense clump of gas in the sur- 



rounding medium. The jet slows, energy is given up to the 
system as heat, and H2 in the impacted gas is sufficiently 
excited that emission takes place in IR vibrational bands. 
C-type shock simulations show that the gas that emits 
most strongly is moving at a velocity of about ten kms~^ 
less than that of the initial jet. For example, there may 
be a jet with an initial velocity of 20 kms^^ impinging on 
gas with a pre-shock density of 10^ cm^'^. In the steady 
state which eventually forms, the heated gas emits most 
brightly in H2 when the flow has slowed to ~10 kms~^, 
according to C-type shock models llWilgenbus et al hood 
iLe Bourlot et al 2002*: iFlower et al.ll2003l L. Kristensen. 
private communication). The gas however continues to 
emit in vibrationally excited states of H2 at even lower 
velocities. 

As the initial collisional process described above takes 
place, more gas is entering the system from the continuous 
jet. This newly arrived portion of the jet encounters gas 
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Figure 12. Left hand panel: a VLT image of flow no. B- 
43, Table EH taken with the NAOS-CONICA adaptive 
optics system, with a spatial resolution of ~35 AU: see 
iLacombe et al.l l)2004l) . Right hand panel: the same ob- 
ject obtained with the GriF instrument on CFHT, with a 
spatial resolution of 70 AU. Contours show velocity and 
colours show brightness in counts per 400 sec. The max- 
imum in velocity is -36 kms^^ at -20.2" E, -6.5" N. The 
interval between contour lines is 4 kms^^. 

which is moving several kms^^ more slowly. Where these 
two portions of the flow meet, there forms a structure 
called the " Mach disk" , which is itself a shock, given that 
the relative flow sp eeds are super- Alfvenic . Thus there are, 
to use the model o effectively two 

shocks, one at the Mach disk, and one around the apex 
of the system where the gas is impacting the ambient ma- 
terial. The gas flows in front of and around the outside 
of the Mach disk, forming a bow structure with so-called 
"working surfaces" at which weaker shocks form in 3D 
around t he region of the Mach disk. This is illustrated in 
Fig. 1 of lRaga fc CabritI |199^. In all parts of this flow, 
the gas cools and emits: hence the structure is seen in H2 
emission. 

This general model flnds support in our observations 
through the fact that the physical apex of the shock, where 
it can be identified, is not where the highest velocity is 
found. The physical tip of the shock has slowed down, 
following momentum conservation taking account of the 
compression of the gas (pivi — ^2^2)- From our observa- 
tions we find that typically the tip of a bow shock is at 
a velocity between 1/2 and 1/3 or lower of the maximum 
velocity. When a magnetic field is present, as here, the 
working surface becomes a greatly extended zone whose 
dimensions may be several hundred AU or more, but 
are typical ly 50-60 AU in high density regions of (say) 
10^ cm-'^ llWikenbus et"ani2000l: iLe Bourlot et al.ll2002t 
iLacombe et al.ll2004|) . 
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Figure 13. GriF data for three bow shocks within the out- 
flow Region B. Left-hand image: H2 emission brightness in 
v=l-0 S(l). Golours denote brightness in counts per 400 
sec. Right hand image: velocity data for the same field. 
The colour scale show velocity in kms~^. The orientation 
of the images is such that north points at an angle of 24° 
measured clockwise as indicated on the figure. The objects 
involved, lying 15" to 20" SW of Source I (see Fig.gJ are 
B-43 to the south (see also Fig. I12|) and B-35-37 in the 
north, see Table IXSl 



Following the above model, the Mach disk is in general 
the point of highest velocity in the shock relative to the 
surrounding gas, assuming a constant jet playing on the 
medium. The question arises, is the Mach disk brighter 
or less bright in H2 emission than the portion of the gas 
in front of it, which is moving more slowly? One might 
suppose that the brightest emission would tend to lie in 
front of the Mach disk, since it is in front of the disk that 
bulk energy is being turned into heat. However in G-type 
shocks, as in OMGl ijVannier et alJ2o'oillKristensen et alJ 
I2OO6I) . the extended shocked zone may well cause the Mach 
disk to merge with the bow region. A quantitative un- 
derstanding of the velocity-brightness structure of G-type 
shocks requires high density models which include both 
chemistry and at least 2D hydrodynamics. At all events 
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no double maxima in the brightness structure of shocks 
are seen in the present observations. 

In the above we assume a stable continuous system 
generated through a stable jet of constant speed. The ex- 
citing jets cannot be readily identified as distinct struc- 
tures in OMCl but where they are observed in class 
and class I object, proto stellar jets may appear time- 
variable in their velocity llBontemps et alJ Il996t IStankd 
l2000t iBallvet al.ll2000l: ' lLarsonll2003|) . This would pro- 
duce a more complex shock structure with several emission 
maxima, but does not in itself involve any shock forma- 
tion mechanism fundamentally different from that out- 
lined above. One point to be noted regarding variable jets 
speeds is that this can result in both forward and back- 
ward facing, or reverse, shocks. Reverse shocks are thought 
to be formed when an initially slow wind becomes a fast 
wind which then overtakes and compresses the material 
ahead of it. This gives rise to clumped material where the 
fast jet impinges on the slower jet. In the limit of a variable 
jet, the jet consists of spatially separated clumps of ma- 
terial each of which individually forms independent shock 
structures. This is illustrated in fig. 1 of .Hartigan (,1,989.) . 

A further point is the relative dimension of the jet 
compared to that of any clump of material on which it 
impinges. In an idealized view, forward bows are formed 
when the jet radius is smaller than the dimension of the 
slow-moving medium and reverse bows form when dimen- 
sion of the obstacle is smaller than that of the jet. 

5.1.2. Determination of orientation angles 

On the basis of the previous section the brightest emis- 
sion and the highest velocity are expected to be closely 
spatially associated. The shock travels in a direction from 
maximum velocity to maximum brightness. This model 
will be referred to as the velocity to brightness model. 
In some cases observations clearly support this model. 
An example is shown in Fig. object B-43, Table 1X31 
Here the maximum velocity is spatially displaced from 
the maximum brightness by 0.3" and the object gives ev- 
ery appearance of a bow shock travelling in the direction 
of maximum velocity to maximum brightness. A counter- 
example to the above is however furnished by object 2-20 
in Table IX!2l Here the bow shape of the emission suggests 
motion in a NW direction, whereas application of the ve- 
locity to brightness model for the shock motion indicates 
motion in the SE direction. Furthermore, in areas with a 
high spatial density of flows, perhaps indicative of a mul- 
tiple star forming region (see sect. 15.2(1 . closely situated 
pairs of outflows where one is red-shifted and the other 
blue-shifted are seen. However, for all pairs of red- and 
blue-shifted shocks where position angles are determined 
by the maximum velocity to maximum brightness flow 
direction model (see below), the flows appear to be con- 
verging in the plane of the sky. If any pair of flows formed 
a bipolar outflow, they would clearly be expected to be 
found diverging, that is, moving apart from one another. 
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Figure 14. An illustration of the proposed geometry by 
which position angles are estimated in the column den- 
sity model (see sect. 15.1."^ . V's marks the position of the 
observed maximum in velocity and B's in brightness in 
the plane of the sky. a a flow moving in the line of sight: 
the maxima in brightness and velocity coincide, b A flow 
moving at an angle to the line of sight: the maximum in 
brightness will appear at some point behind the tip of the 
shock. A vector through the two maxima then gives the 
direction of the flow projected onto the plane of the sky, 
the position angle. 



An example of this discrepant behaviour may be found in 
Fig. ITSk showing flows 1-76 and 1-81 (see Table rO|l . 

Failure of the velocity to brightness model will arise if 
effects of optical depth are paramount in determining the 
location of maximum brightness. The brightness of emis- 
sion is a function not only of the local efficiency of excita- 
tion associated with local loss of kinetic energy, but also of 
the column density of gas which is observed. Thus a shock 
travelling at some arbitrary angle may appear brighter 
behind the position of maximum velocity by virtue of its 
bow form, although locally the emission is weaker than 
at the tip of the shock, where the latter is of low column 
density. This model is illustrated in Fig. ^1 and was de- 
scribed in Paper I. This model will be referred to as the 
column density model. 

The common feature of both models is that the max- 
imum in emission and the maximum in velocity lie on a 
vector defining the orientation of the shock in the plane 
of the sky. However, in the absence of a full 2D or 3D 
C-type shock code for H2 excitation including chemistry, 
it is not possible to determine which processes and thus 
which model may be relevant in any particular case, save 
where the morphology clearly indicates a bow shock (as in 
Fig.lT^. Tables IXtI - 1 A . 31 a note position angles measured 
east of north (i.e. counter-clockwise) based on the veloc- 
ity to brightness model where shocks are travelling in the 
direction from the position of maximum velocity to the 
position of maximum brightness in the plane of the sky. 
However, with a few exceptions where the morphology of 
the ffow would appear to dictate the absolute position an- 
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Figure 15. Number of flows in each of the three main regions with position angles in a given range (see also star 
diagrams in Fig. 0J. 



gle, e.g. fig. El these position angles are uncertain by 180° 
with respect to absolute direction as discussed above. The 
exceptions are marked with a star and are found only in 
Table I A. 31 for the outflow region B. The uncertainty in 
the position angle is estimated to be ~ ±25°, decreasing 
with increasing displacement. The uncertainty stems from 
the determination of the exact positions of the maxima. A 
displacement of <0.25" of maximum brightness and max- 
imum velocity is close to the resolution limit, and small 
variations in displacement would have a large influence on 
derived position angles. Position angles are therefore only 
determined where the displacement is >0.25". 

An important consideration is whether the position an- 
gles are dominated by a general outflow motion from the 
BN-IRc2 area or are approximately isotropic in distribu- 
tion in each main region (Peak 1, Peak 2 and Region B). 
Fig. El show the distribution of position angles binned 



over 30° intervals, with the reservation that individual po- 
sition angles may be incorrect by 180° (see above discus- 
sion). In this connection, the assumption has been made 
that one or the other model for position angles (veloc- 
ity to brightness model or column density model) will 
dominate in a given zone, given that a single mechanism 
causing shocks predominates in any one zone, that is, 
outflows from star forming regions or a large scale flow 
impinging on the region (as in Region B). With this as- 
sumption, any preferred directions can be found by look- 
ing for anisotropics in the distribution of position angles. 
Rebinning the histograms in Fig. 1151 to just two bins 
of width 180° and then shifting the absolute positions 
of the bins to cover all possible combinations will show 
anisotropics as large differences in the number of flows in 
each bin. This corresponds to placing a line across the star 
diagrams in the inset in Fig. 0] and counting the number 
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of flows on each side of the line. The greatest difference 
between the numbers above and below the line will ap- 
pear when the line is perpendicular to the direction of 
an overall motion, while there will be roughly equal num- 
bers when the line is parallel to any overall motion. If the 
line is rotated (corresponding to the bins in Fig. 1151 being 
shifted) any preferred direction will appear as a maximum 
in the ratio between the number of flows on each side of 
the line at an angle perpendicular to the outflow direc- 
tion. If there is no preferred direction of motion the ratio 
should be approximately unity independent of the angle 
of rotation. Fig. El shows this ratio as a function of an- 
gle measured counter-clockwise from north for each main 
region. The ratio is calculated so that it is always >1. 

Region B (Fig. 116b ) shows a high degree of anisotropy 
around 140-160° and a minimum around 50-70°. As de- 
scribed above, this property with a maximum and a min- 
imum perpendicular to the maximum, indicates that the 
motions are dominated by an overall outflow motion. In 
this case the motion is along a line oriented ~ 50-70° 
counter-clockwise from north, that is, approximately NE- 
SW. The nature of this outflow in Region B will be dis- 
cussed further in sect. El 

Fig. 116b shows a minimum around 130-160° and a 
maximum at 20-50°. This is similar to the structure ob- 
served for Region B although the maximum is less pro- 
nounced. This suggests that the observed flows in Peak 1 
are also partly described by an outflow motion (towards 
the NW in this case), but not to the same degree as in 
Region B. Finally Peak 2 fFig.ll6b) shows only very weak 
maximum-minimum structure, with a distribution of po- 
sition angles which is close to isotropic. This suggests that 
the flows observed in Peak 2 are only weakly affected by 
general outflow motion and tend rather to be individual 
fLows, such as might arise from local star formation. 

5.2. Protostars and Star forming regions in OMCl 

Patterns of gas motion and spatial distribution of flows are 
examined here for Peaks 1 and 2. Our discussion provides 
further evidence for on-going low mass star formation in 
OMCl. The major source of energy injection into OMCl 
is at 0.1 pc (20000 AU) scales through outflows from mas- 
sive stars an d as part of a large sca le cascade of turbu- 
lent energy ijCustafsson et al.ll2006tl . Thus many of the 
flows, both large and small scale, identified in the present 
work are an expression of the turbulence inherent in the 
gas thro ugh massive stellar outflows from the BN- IRc2 
complex ( Shuning et a l. '2004; Crccnhih et al."2 0o'4bl) and 



from source I (see below: iGreenhill et al. 2004ai1 . However 
the gas dynamics may also be influenced by the events 
which we focus upon here, namely low mass star forma- 
tion. This would provide an additional source of energy 
injection and generates flows at smaller scales of a few 
arcseconds, 1000-2000 AU, through bipolar outflows from 
isolated protostars or multipolar outflows from multiple 
star formation (see Figs.lTHk-c'). 



Table 1. Good protostellar candidates, "candidates", 
(upper half of table) and weaker candidates, "possible 
candidates" , (lower half) in the area of O MCl covered in 
the pr esent work, based on data taken from lMuench et al] 
1)200^ . For purp oses of cross-re f erence , we have adopted 
the numbering of iMiiench et al.l 1 2002|) . 

Protostellar candidates in OMCl 



K-L > 1.5 



Candidates 
J-H > 1.0 or J- 



H undetermined. 



Quench ct al. 


Coordinates 


Colour indices 


numeration 


rel. to TCC0016 


K-L 


J-H 


00506 


-18.3 


-17.1 


5.77 




00543 


9.5 


-7.2 


3.29 




00557 


-0.9 


-5.0 


2.32 


>1.82'' 


00560 


-38.4 


-2.3 


3.03 


>0.82'' 


00569 


-24.0 


-0.1 


3.57 




00572A 


-3.6 


0.6 


1.52 


>2.67'' 


00576 


19.2 


1.6 


1.77 


1.07 


00578 


-14.6 


1.6 


4.42 




00580 


4.4 


2.4 


2.76 


3.77 


00583 


-8.3 


3.0 


1.50 


1.86 


00596 


14.6 


5.9 


1.53 


1.11 


00598 


-8.4 


6.3 


2.99 


3.30 


00603 


-0.6 


7.5 


2.22 


1.02 


00614A 


-2.9 


9.4 


2.81 


>4.65" 


00645 


-17.6 


17.1 


1.82 


1.32 


00709A 


-9.2 


34.7 


1.99 


1.25 


00726 


-16.7 


39.4 


2.55 


2.94 


Possible candidates 



No L band ; H— K > 2.0 or K band photometry only. 



Jvlucnch ct ^ Coordinates 
numeration rel. to TCC0016 


H-K 


M 


K 


00552 


-20.9 


-6.1 


2.81 


13 


03 


00554 


-4.7 


-5.8 


>2.79'' 


15 


01 


00561 


15.8 


-2.1 


3.43 


13 


22 


00570 


-5.7 


0.0 


>3.00'' 


14 


80 


00585 


-10.7 


3.4 


3.19 


10 


96 


00591 


-11.1 


4.9 


3.49 


10 


97 


00594 


1.5 


5.2 


>3.85* 


13 


95 


00609 


-40.8 


8.3 


>1.63* 


16 


17 


00610 


-9.2 


8.5 


6.09 


12 


44 


00634 


-8.6 


14.2 


3.35 


10 


78 


00665 


-17.1 


21.7 


2.04 


11 


01 


00691 


-25.1 


29.3 


3.05 


13 


39 


00695 


-37.4 


30.4 


3.13 


15 


59 


00698 


-44.0 


31.1 


>2.9l'' 


14 


89 


00713 


-47.7 


36.1 


>2.99* 


14 


81 



Muench et al 



quote a detection limit of 18.15 in J band, 
quote a detection limit of 17.8 in H band. 



We first consider the IR observational evidence for 
the presence of early protostars in OMCl. To this end 
we have compiled a table of protostellar candidates 
in OMCl (table HI). We u se the criterion reported in 
iKenvon fc HartmanrJ (^^B); for YSOs in Taurus, that an 
object with K— L > 1.5 is almost always of a protostellar 
nature (specifically class I) and assume that this also holds 
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Figure 16. The ratio between the number of flows on each side of a hnc through the star diagrams in Fig.Qlas a function 
of the angle between the Hne and north measured anti-clockwise from north. A high value suggests an anisotropy in 
the distribution of position angles across the line while a value of unity indicates an isotropic distribution. 



fo r Orion. W e apply this crit erion to the photo metry data 
of iMuench et al. ( 20q3) and lLada et~an pOO^) for Orion. 
We furthermore use the criterion that J— H m ust be >1.0 
to exclude, for example, unresolved binaries ijLada et alJ 
[2000'! . The resulting protostellar candidates are shown in 
tabled For e ase of reference we hav e adopted the number- 
ing system of lMuench et all l)2002|) . Where a ">" symbol 
is encountered, this implies that J band data are absent, 
in the upper section of the table, or H-band are absent in 
the lower. 

We have divided objects into those which are very 
likely protostars, labelled "candidates" (upper part of ta- 
bleOJ and those where the identification is more tentative, 
labelled " possible candidates" (lower half of table . In 
addition there are objects for which the colours suggest 
that they may be excluded as protostars. These are not 
shown in table This last category forms the majority 
of point-like objects within our complete field. However 
within Peaks 1 and 2, these objects occur in roughly 
equal numbers to the sum of "candidates" and "possi- 
ble candidates". Fig. 1171 shows the positions of all objects 
in the 3 categories. Candidates and possible candidates 
have been marked with circles and triangles respectively 
with their accompanying number for reference to tabled 
Non- candidates have been marked with a squar e. From 
a comparison with fig. 1 in IShu oing et al.l l|2004l) we find 
that among the objects in tableware source n (00598), 
IRc4 (00585), IRc7 (00610), IRcll (00614) IRcl2 (00603) 
and IRcl5 (00634). Furthermore IRc9 is found to be iden- 
tical with 00726. Source n is discussed in more detail in 
sectional IRc9 will be mentioned in the following section, 
sect. 15.31 

The term "candidates" covers objects which have K— L 
> 1.5 and which have either J— H > 1.0 or where J— H is 
undetermined. The term "possible candidates" covers ob- 
jects with H— K > 2.0. As far as the authors are aware, the 
latter is a new criterion for the identification of protostars, 



motivated by a lack of L-band photometry. The new cri- 
terion has been established as follows. An extensive anal- 
ysis of the H-K colour index of the stars in lMuench et alJ 
1)20021) has been carried out. It has been found that for 
objects with K-L > 1.5, 39% have H-K > 1.5 and 14% 
have H— K > 2. By contrast, for objects with K— L < 
1.5 - which are not protostellar candidates - less than 
2% have H— K > 1.5 and only one object (0.3%) has 
H— K > 2. On this basis a H— K colour index of > 2 is 
adopted as a tentative identification of "possible proto- 
stars" in table ^ and in Fig. El We note that the value 
of H— K cannot however be used to exclude protostel- 
lar candidates. In this sense the objects listed in tabled 
form a minimum set. Under " possible cand i dates" in ta- 
ble ^ we also list object from lMiiench et~ ] 1I2OO2I) which 
have only been detected in K band. These are shown 
with H— K> some value. The K-band detection, without 
J and H, suggests that these objects are deeply buried in 
the dusty gas. We conclude that there 17 good protostel- 
lar candidates in our field with a further 15 (or more) 
possible candidates. A comparison has also been made 
with the Chandra Orion Ultradeep Project (COUP) X-ray 
data (Getman et al. 2005; Grosso et al. 2005) and sources 
in tabled COUP contains numerous point sources in the 
direction of the region studied here and a detailed compar- 
ison will be made elsewhere. Suffice it say at present that 
there is a remarkable correspondence between candidates 
in table ^ with COUP sources with 14 out of 17 spatial 
overlaps, whilst the so-called possible candidates are anti- 
correlated, wi th only 3 out o f 15 sp atial overlaps. Overall 
about half of lMuench et al ] l|2n02l) sources are spatially 
coincident with COUP sources. 

We now turn to the positional correlation of these ob- 
jects with possible outfiows and H2 emission in general. In 
some regions two closely spaced shocks moving in opposite 
directions are seen, which could indicate that a protostar 
is situated in between, causing the shocks. An example 
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Figure 17. The positions of point-like sources from iMiiench et aD l(2002l) in the observed field of OMCl. Objects 
have been divided into 3 categories based whether they are candidates (circles), possible candidates (triangles) or 
non-candidates (squares). Numbers refer to the labels in table 



of this is shown in Fig. llSh . where a pair of flows are 
found, one rcd-shiftcd and the other blue-shifted. Data 
in table IXTI show that position angles of these flows are 
roughly 180° apart, noting the indeterminacy with regard 
to their relative direction. However, if the positions of pro- 
tostellar candidates in tableware compared with the po- 
sitions of candidate bipolar outflows, there is no clear cor- 
relation between the two. In only 2 positions, 15.75" E 
-2.1" N (object 00561, see table HJ and 1.50" E 5.2" N 
(00594) a candidate protostar appears directly over a pair 
of flows, whereas in many cases a "bipolar pair" is found 
where there is no evidence of a protostar. Thus a "bipolar 
pair" in itself does not appear to be evidence of the pres- 
ence of a protostar. However the caution should be made 
that we are considering the earliest stages of protostellar 
development (class O/I). This phase is characterized by a 
very high degree of obscurat ion, an d new, deeper obser- 
vations than those of Mucnc h"et al.l l)2002l) . used here to 
detect the protostars, may yield different results. 

As described in sect. 14.41 there is a marked difference in 
how flows are distributed in Peak 1 and Peak 2. In Peak 1 
the flows are mostly uniformly distributed across the area. 



whilst in Peak 2 the shocks often appear in small groups of 
6-8 shocks within an area of ~5" in diameter (see Fig. ITTI 
and Figs. [T5b and c), which corresponds to ~2400 AU, 
at the distance of Orion. Since most stars are expected 
to form as part of a binary or multiple system ijPudrit d 
2002; Larson 2003,) , it seems plausible that these groups 
of shocks indicate places where multiple star formation 
is ongoing: the presence of 2-3 protostars would account 
for 4-6 shocks from bipolar outflows. The flows in these 
regions are not however generally found as pairs of shocks 
moving in opposite directions, as may be seen in Figs. 118b 
and c with reference to the position angles in Table IX!^ 
This supports a picture in which multiple star formation 
is a complex process involving interactions between the 
protost ars and disks causing jets which are episodic or 
pulsed l)Larsonll2003l and references therein). 

We now examine whether we can associate regions with 
large densities of flows, as just mentioned, with protostars 
listed in tabled Moving from east to west (left to right) in 
fig. llll the first region of high flow density is around coor- 
dinates 7" E 0" N. Only one stellar-like object is observed 
along the line-of-sight to this region: the proplyd 154-240 
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Figure 18. Various areas with multiple flows in a small area: a A possible bipolar outflow. Flows no. 1-76 and 1-81: 
see Table lA.ll b A complex flow pattern. Flows no. 2-42, 45, 46, 47, 48, 49: see Table IA.2| c Another example of a 
complex flow pattern. Flows no. 2-16, 18, 19, 20, 24, 25: see Table fST2l Colours denote the brightness in counts per 
400 s exposure. 



l)0'Dell et alJll997HSchultz et al.lll999(l . Since this object 
has been observed in the visible by HST it cannot be a 
protostar buried in the dense, dusty gas. 

Moving away from a line-of-sight positional correlation 
of possible protostars and H2 emission, there is a good 
protosteUar candidate ~4" NW or 1600-2000 AU in the 
plane of the sky. This is the object 00580 ^ which has 
K— L, J— H > 2 fulfilling our criteria for a buried protostar. 
This is therefore a good candidate to power at least some 
of the emission that is observed in this highly active zone. 

The second region of high flow activity in fig. ^2 is 
around coordinates 0" E 6" N. Within a radius of a few 
arcseconds there is one star which satisfies the protostel- 



lar criteria, namely 00603 with K-L = 2.22 and J-H — 
1.02. There are two other stars for which we do not have 
L band data available and therefore we are unsure of their 
nature. One of these is 00594 which is listed in tablenun- 
der " possib l e cand idates" This object li e s just outside the 
iLada et all l)2004|) frame. iMuench et alJ l)2002(l do not de- 
tect this object in J and H and record niK = 13.95. Thus 
it is a highly obscured object and therefore accordingly 
appears in tabled A further object is 00606. This object 
has no L-band data, and with H— K = 0.57 this is not 
an obvious protosteUar candidate. However, one should 
recall from the above discussion that the H-K colour in- 
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dex is useful only for including, not excluding, protostellar 
candidates. 

We turn now to the third flow concentration, which 
is centred on —22" E 37" N. Surveying the region 
within a radius of ^ 4", we find the star IRc9 
(00726) fe mith Ra,ihll2nn,'j) . 



Ballvl l20n,'iD . with J-H and K-L both 
>2, fulfilling the protstellar criteria, see tabled There is 
also another object (00715) within a few arcseconds from 
this concentration of flows, but again there is insufficient 
photometry to make an assignment as to the nature of 
this object. 

Looking at Peak 1 in general, rather than at concen- 
trations of flows, there are a further two protostellar can- 
didates, 00645 and 00709A, both situated along the edge 
of the H2 emission seen in Peak 1. There are also another 
5 possible protostellar candidates (00665, 00691, 00695, 
00698, 00713) which are spread across Peak 1. In Peak 
2, there are a further 3 protostellar candidates (00543, 
00576, 00596), and one possible candidate, 00561, associ- 
ated with the very bright H2 emission at 16" E -2" N. The 
object 00561 has in fact an associated 140 AU diameter 
object observed in the L-band with the VLT (Paper I), as 
mentioned in sect. 13.21 We note that with the exception 
of 00561 (in Peak 2) no stars or protostars in Peak 1 or 
2 are directly associated with the brightest areas of H2 
emission. 

In region B there are 8 protostellar candidates and a 
further 8 possible candidates. Half of these are situated in 
the dense core of OMCl around Source I and may con- 
tribute to H2 emission in this region. The other half are 
widely distributed over region B. 

The discussion here and in preceding sections points 
towards Peak 2 being a more active star forming region 
than Peak 1 with similar numbers of protostellar candi- 
dates despite Peak 2 being only half the size of Peak 1 in 
the plane of the sky. This is an unexpected result given 
the proximity of the two regions. One possible explana- 
tion is that the ionization front from O^Oii C has cleared 
away most of the diffuse gas in Peak 2, leaving only the 
densest parts of a previously more extensive cloud and al- 
lowing us to see the dense star forming cores which would 
otherwise have been obscured fro m view. T his interpreta- 
tion is supported by the work of IPoeee et al. (1992) and 
iTakami et al.l ^^2002^ who show evidence of an ionization 
front just north of Peak 2. 

The above considerations ignored the presence of 
WRVs in Peaks 1 and 2. In Peak 2 WRVs form only 20% 
of the observed flows but in Peak 1, they form more than 
50% of the flows. However analysis presented in sect. 14.31 
showing the organised nature of WRV flows in Peak 1, 
suggests that WRVs form a subset of flows which should 
not be associated with local star formation. This is further 
discussed in the section directly below, sect. 15.31 

At the start of this section we asked whether H2 emis- 
sion might in part be caused by the presence of protostars 
with associated outflows. Our data show that this is the 
case but we do not demonstrate that the majority of bright 
emission is due to protostellar activity in OMCl. By their 



nature the most important sources may be deeply buried. 
However for the present we conclude that the emission 
arises from a mix of protostellar and massive star activity. 
In particular in Peak 1, present data suggest that emission 
may be dominated by massive outflow from the vicinity 
of Source I and the BN region. 



5.3. The origin of weak radial velocity structure in 
Peak 1 

The spatial distribution of weak radial velocity compo- 
nents (WRVs) in Peak 1 revealed a band, ~2-3" in width 
and more than 30" long passing diagonally across emis- 
sion in Peak 1 (see sect. 14. 3|) . It was also noted that 2/3 
of these WRVs showed a much more broken-up structure 
than other features, including other WRVs. This section 
briefly explores possible scenarios for the origin of this 
band of WRVs. 

The first scenario is that the band of WRVs may rep- 
resent a less energetic part of the outburst which created 
the "bullets" or "fingers" l)Allen fc Burt on"l993'). These 
are clumps of gas moving radially outwards at high veloc- 
ities of up to several hundred kms~^ ijPoi et alJl2002l|) and 
appear to hav e an ori gin in the area around BN-IRc2 an d 
Source I ( Sa,l a,s et abl ll 999: O'DcU 2001; Do i et a]Jl20ol . 
If the WRVs are moving at a velocity of 40 kms~^ - the 
upper limit of motions observed here for other clumps of 
H2 - in 1000 years they would be found at a position 
about 18" from the BN-IRc2 complex, essentially where 
they are observed. The typical mass of each WRV can 
be estimated to be ~ 1.5 x 1O~^M0, given a density of 
10^ cm^'^ and a diameter of 500 AU. This yields a total 
bulk energy for outfiowing hot H2 in the WRVs of the or- 
der of 5 X lO^'' ergs, for 20 such objects. This is about 
1% of the energy ass ociated with th e bullets whose prop- 
erties are reported in lBurtonll|l997l) . 



A second scenario is that the band of WRVs is the re- 
sult of bright photodissociation zones powered by O^Ori C, 
but with enhanced brightness due to geometrical effects 
(see below) . The morphol o gy of the ionization fr ont, de- 
scribe d in IWen fc O'Delll ^ll995^ and refined in lO'DeJ 
shows that the bright H2 clumps observed in the 
present work arise fr om dense gas within t he Hii region 
created by 9^0ri C ijKristensen et al.1l2003l) . Ionizing ra- 
diation has cleared away the diffuse gas and is now slowly 
eating its way through the denser clumps of gas remain- 
ing. Thus the line of WRVs along the eastern edge of Peak 
1 may be the relic of a denser region that existed in the 
original molecular cloud. Whilst the brightness of the H2 
emission cannot be reconciled with a face-on PDR (see 
sect. 13.2(1 ■ PDR emission may be en hanced local l y thro ugh 
a geometry originally proposed in iField et alJ ^^i)- In 
this model, a choppy surface of H2 yields chance lines-of- 
sight with high column density. This is consistent with the 
broken up appearance of 2/3 of the features in this zone, 
mentioned in sect. 14.31 



22 



H. D. Nissen et al.: Shocks in Orion 



If we consider only those WRVs lying in the band of 
width 2.5" between points P and Q in fig. ^ a further 
possibility is that these may be associated with an out- 
flow in the plane of the sky from IRc9 (00726 in table Q 
and fi g. [TTIl. IRc9 is l ocated ~2" from the centre of this 
band. lSmith fc Ballvl l| 2005>> show that IRc9 is a YSO pos- 
sessing a clearly defined disk. The major axis of the disk 
is oriented at right-angles to the WRV band which sug- 
gests that WRVs in the band could represent a highly 
collimated outflow in the plane of the sky from IRc9. In 
sect. 14.31 we found that there was a significant overrepre- 
sentation of WRVs in Peak 1. If the analysis is repeated 
but without including the WRVs in the band, a significant 
overrepresentation of WRVs in Peak 1 is still found: based 
on the same model used in sect. 14.31 the expected number 
of WRVs in Peak 1 is ~12 while 32 WRVs are actually 
observed. Thus a significant discrepancy remains. 

6. The outflow zone, Region B 

6.1. The geometry of the region: radio observations 
and the present data 

An outflow in Region B has previously been identified in 
radio-observations involving both proper motions and ra- 
dial velocity measurements of H2O, OH and SiO masers 
and thermal emission in SiO (v=0) and NH3. Evidence 
of shocked gas within the flow is provided through ob- 
servations of emission from high rotati onal states of SO, 
SO2, H ON, where data are reviewed in lGenzel fc Stutzk] 
l)l989tl . The outflow has also been identified in sub sequent 
observations reported in lWright et al ] |ll992Lll99fil) in data 
involving 16 molecu lar species, in particular SO and SO2 
l|Wright et al.ll996|) . The outflow is associated with source 
I whose location is given in Fig. [T^ 

The expansion velocity associated with the above- 
mentioned outflow in Region B on the basi s of radio ob- 
serva tions was fou nd to be 18d=2 kin s~^ l|Genzel et all 
Il98ll) . The data of iGenzel et"al] l)l98l[l showed a flow ex- 
tending ~20" SW of the BN-IRc2 region, which is very 
simil ar to the SW exte nt of the H2 emission recorded 
here. lGenzel et all lll98l|) also note, through their observa- 
tions of numerous objects at closely similar velocity, that 
the flow appears to have been stable for >3000 years. 
This would seem to preclude the scenarios discussed in 
IShuping et alJ l)2004l) in which source I may be the origin 
of both of the major outflows in this region, through un- 
stable behaviour, such as flipping of the orientation of the 
ou tflow. 

iMenten fc ReidI l)l995|) performed SiO maser observa- 
tions of the centre of the region around source I, extend- 
ing over ~50AU (0.11"), and found both blue- and red- 
shifted flows, with the blue-shifted component showing a 
visr between —13.0 and —6.2 kms~^ corresponding to a 
blue-shifted flow relative to the ambient gas of between 25 
and 18 ±6 kms~^, taking into account the uncertainty in 
the visr- The centre of the SiO maser emission was found 
to be located at the position of source I. The same spa- 




10 -12 



Figure 19. The positions of BN, radio source I, radio 
source n and IRc2 a-d, superimposed on velocity inte- 
grated emission in the v=l-0 II2 S(l) line. 



tial relationship of SiO masers to source I is recorded in 
iDoeleman et al ] (IT999). Further recent SiO observations 
covered an inner area of <2" diameter (900 AU) and in- 
volved both SiO maser and non-maser emission, the latter 
in v=0, and maser emission in H2O , inclu ding proper mo- 
tion studies ijGreenhill et alJ20o'4aj) . As in lMenten fc ReidI 
|1995.) . the expansion centre of these motions was found 
to be coincident with source I, in this case to within 30 
mas (^15 AU). Radio source n (hereafter simply source n, 
see Fig. I19II a few arcseconds fr om source I is also s ignif- 
icant. From the obser vations of iGenzel et al ] lll98lh and 
iMenten fc ReidI (|l995') source n is found at the centre of 
expansion of a shell of H2O n iasers, perhaps const ituting 
an independent outflow zone l|Shuping et alJl2004j) . 

Our present observations measure gas motions relative 
to the ambient surrounding gas and therefore may also be 
interpreted as a radial component of an expansion velocity. 
The expansion velocity of 18i b 2 km s~^ derived from the 
observations of IGenzel et alJ lll98j) matches accurately 
the mean velocity of —18 kms~^ observed here for the 
hot H2 in Region B: see Fig.jS] The veloci ties measured in 
this w ork a re also consistent wi th data in lMenten fc R.eidI 
lligg.-j) a,nd lWnght ellil l|l996h . 

The much debated issue of the origin of large scale 
ou tflows in OMC l , exa mined in detail for example 
in IMenten fc ReidI l)l995|) and elsewhere, has recently 
achieved a new perspective wi th high spatia l resol ution 
(0.3-0.5") IR ineasurem ents in lShuping et alJ l)2004|l and 
iGreenhill et al] ll2004b() . These observations show that 
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there are a number of potential YSO candidates buried 
within a zone in the north-east of Region B. IRc2 was long 
believed to be the main power source in OMCl and the 
source of the outflow activity in the regio n. However these 
and other high resolution observations ijPougados et al.l 
have resolved IRc2 into several individual so urces. 
Thes e may or may not be self-luminous (Shuoing et al.l 
|2004}| . At all events the current most likely candidates for 
outflow sources are sources I and n. 

The positions of sources are shown in Fig. [T^ 
Sour ce I (-5.95" E 8.56" N), is a deeplv buri ed O- 
star l|Menten fc Reidlll995t ICreenhifl et al] l2004a»). It is 
the most luminous source in the h ot core of OMCl with 
a luminosity of 5x1 0^ to IO^Lq l)Menten fc Reidlll995t 
iBeuther et al.l lioO^ and see sect. Iti.2|l . c ompared for 
exam ple to ~25OOL0 for the BN-obiec t llGezari et al " 



or ^2000 L0 estimated for source n ijCreenhill et al 
2004bl ). Source I is obscured in the infrared and cannot be 



detected at wavele ngths as long as 22 um { ^rccnhill ct alj 
l2004bj) . Source n jcreenhill et all 



2004b 1 is located at 



Shuping et aD l|20a 



-8.35" E, 6.40" N and according to 
may give rise to an outflow in this region (see above and 
sect. l6.3ll . 

Base d on their data for SiO a nd H2O emission around 
source I. iGreenhiU et all l|2004aj) proposed a model with 
an accretion disk and a wide angle outflow pointing NE- 
SW wit h an opening an g|le of ~90°. This model replaces 
that of iGreenhill et alJ l[l998'). The morpho logy of the 
outflow from source I, newly proposed in lOreenhill et all 
l)2004a|^ . is shown in Fig. 1201 lower left panel. 

The axis of the outflow from source I passes through 
a region where we find a striking II2 emission feature. 
Velocity data for this zone is shown in the upper right 
hand frame of Fig. EH Three fast jets (flows no. B-7, 
B-9, B-10: see Table IA.3|l are seen impinging on the 
surrounding medium. A VLT image of the same area 
ijLacombe et al.l f2004l) is shown in the lower right-hand 
frame of Fig. 1201 The model of the outflow from source I 
is superimposed upon it, to show the spatial relationship, 
where the accuracy of the relative positions is ~0.2". The 
position of the maximum brightness of this outflow re- 
gion, measured in H2, is 1.84" W and 0.96" S of source I. 
The appearance of the H2 emission in this region suggests 
that this is a zone where the radio-detected outflow bursts 
through the dense gas surrounding source I. Thus our ob- 
servations show the IR counterpart of the 0-star outflow, 
hitherto only obs erved at radio-w" a yelengt hs. 

The model of lOreenhill et all l)2004a|) of the outflow 
from Source I is of an 0-star with an edge-on disk and 
a wide-angle (90°) bipolar outflow in the plane of the 
sky. This begs the question of why no II2 emission is 
observed as part of the outflow to the NE of source I. 
We would suggest that the flow axis may in fact tip out 
of the plane of the sky by ^45°. Flows would then tend 
to be rendered invisible tho ugh dust obscur a tion in this 
very dense zone. Moreover, lOreenhiH et all j2004ah ob- 
serve a substantial overlap of H2O and SiO masers along 
the line-of-sight. However these two molecular masers 



trace different gas den sities and temperatures, and as 
ICreenhill et alJ l|2004al) note should not be found within 
the same mass of gas. H2O masers may have associated 
kinetic temperature of 200 to 2000 K, with a preference 
for temperatures around 600 K, and densities of no more 
than a few x 10^ cm~^ but most likely with densities 
<10^ cm"3 l|Yates et al.lll997l) . SiO masers require SiO 
to remain in the gas phase, with associa ted temperatures 
of >1500 K and densiti es of 10^° cm'^ ijPoel et al.lll995t 
Huniphrc vs et alJll996|) . If, however, the outflow cone has 
a significant component in the line-of-sight, as we suggest 
here, this discrepancy may be resolved: the H2O masers lie 
further from the 0-star than the SiO masers, but overlap 
in the line of sight. Further, blue-shifted flows are observed 
to the north of Source I. This is consistent with an outflow 
originating from source I, but tipped out of the plane as 
we suggest. This is illustrated in Fig. |^ - described in 
more detail below - where any material to the left of the 
line-of-sight through source I would effectively appear as 
blue-shifted material to the north of source I. 

In conclusion wc agree with lOreenhill et al.l ( 2004ah 
that the direction of the outflow in the plane of the sky 
is NE-SW. Following the result of sect. 15.11 the position 
angle of the axis of the cone of the outflow from source I 
is '^240°±10°. We suggest however that the flow, rather 
than lying in the plane of the sky, lies at an angle of as 
much as 45° to this plane. 

A schematic diagram of the region, not to scale, based 
upon this model, may be found in Fig. [21 This shows the 
region as if observed from right-angles to the line-of-sight, 
from a position in Peak 2, looking west. Starting from 
source I at the base of the diagram, first SiO and then 
H2O masing regions, overlapping in the line-of-sight, are 
found in the dense gas close to source I. Moving further 
from the 0-star, the outflow breaks through the dense gas 
and dust surrounding source I, generating H2 emission 
features such as that shown in Fig. [5D| and, moving still 
further, other blue-shifted features shown in Fig. 0] and 
listed in Table EH 

The model proposed above requires that the directions 
of flows of H2 emitting regions should show asymmetry in 
the NE-SW direction. Referring to Fig. E| and Fig. El the 
star diagram for Region B in the upper left panel shows 
a 3:1 anisotropy about a dividing line running roughly 
perpendicular to the outflow direction. In this connection, 
recollect that the absolute direction of the flow, ±180° 
is not determined, as noted earlier. The only additional 
assumption at this juncture is that one or other model on 
which the flow direction is based is dominant, see sect. 15. II 



6.2. Physical properties of the outflow source 

This section considers some of the governing parameters 
of the 0-star as s ociated with source I and with region B. 
ICreenhill et all l(2004lJ ) found that source I could not 
be detected at 22 /xm. They concluded that the optical 
depth at 22 nm must be >300. Using the relationship be- 
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Figure 20. Source I. Top left: Region B recorded in the v=l-0 S(l) line of H2. The yellow star shows the position of 
BN and the red star shows the position of source I. The red circle is source n, and the yellow circle is IRc2 d (compare 
to Fig. I19|l . The black outline shows the position of the H2 emission fe ature shown in the right hand panels. Bottom 
left: the geom etry of the wide-angle outflow from source I taken from iGreenhill et al ] l)2004a|) . Bottom right: VLT 
NACO image ijLacombe et al.lE004j) of the outflow region (the "fireball") just SW of source I, outhned in black in 
the top left panel. The inset shows that the outflow from source I (not to scale) points directly at this region. Top 
right: velocity structure of the fireball, showing three fast jets impinging on the medium (flows no. B-7, B-9, B-10: 
see Table IX3|) . 



tween extinction and wavelength given in lRosenthal et aH 
l)2000|) . the NIK optical depth will then be >570. This im- 
plies that the column density of H2 is > 1.5 x 10^^ cm~^. 
This is com parable with the figu re of <8.5 x 10^"' cm^^ 
estimated in lBeuther et al.l l)2004(l . Depending on the av- 
erage gas density around the obscured 0-star, this object 
must be buried several arcseconds within OMCl. For ex- 
ample if the average density is 10^ cm~^, the 0-star is 
buried to a depth of 2" and the mass of gas surround- 
ing the 0-star is 9 M©. This is comparable with the 
estimate of the mass of the 0-star itself of <10 M© in 



Beuther et al. 



Shuping et al 




or the figure of 10-20 M© given in 
These masses are consistent with 
a luminos ity of a few t housand to ^ 10'* L©, consistent 
also with iMenten fc Re id (1995). The most recent esti- 
mate places the luminosit y of source I at ^5x10^ L© 
l|Beuther et al.ll2006ll200^ . 

The mass outflow rate in Region B may be 
timated as follows. Typical gas densities in the 
clumps in Regio n B ha ve been investigated in detail in 
iKristensen et al.l l)2006f) and post-shock values range be- 
tween 10^ to 10^ cm~^, equivalent to an average density 
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Figure 21. The outflow from source I, showmg the appearance of the outflow and its relationship to other objects in 
Region B, seen at right-angles to the line-of-sight from a vantage point in Peak 2 looking towards the west. 



of '-^lO^^^ gcm^'^. Assuming a continuous outflow and a 
measured velocity of 18 kms~^, this represents a mass flux 
of 2 X 10~^^ gcm~^ s~^. The total area in the plane of the 
sky covered by the excited H2 may estimated by using a to- 
tal area for each clump given by the area around each flow 
where the emission is >50% of the maximum brightness 
of the flow. The total area is found to be ^^5 x 10'^^ cm^, 
giving a mass flow of ^1.5 x 10"'^ Mq yr~^. Including the 
red-shifted part of a bipolar outflow, this gives a total out- 
flow rate of ^3 x 10^'^ M0 yr^^. Since only excited gas has 
been included this is a lower limit. This figure for t he out- 
flow r ate is comparable to estimates given in Sheoherdl 
l)200fj) for stars of the luminosity associated with source 
I. 



An independent estimate of the outflow rate may also 
be performed. If the luminosity of the source is accretion 
luminosity and the surface temperature is (say) 1500 K, 
typical of a protostar, then the stellar radius is '^5 AU, 
given the figure of 5x10"^ Lq. The mass infall rate is then 
0.02 Mq yr~^. For low mass st ars the outflow rate is gen- 
erally -10% of the infaU rate l)Richer et al.ll2000() . If the 
same proportion is adopted for the case of this massive 
star, the result is a mass outflow rate which is essen- 
tially the same as that which has been estimated from 
our observations. The corresponding outflow momentum 
rate is ~0.05 MQyr"-'^ kms^^ and mechanical luminosity 
is ~80 L0 . These yalues are tvpica l of a late- type 0-star or 
early- type B-star ( SheDherdl200f)l) . Together these charac- 
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teristic properties combine to make a consistent picture of 
an outflow from source I where this source is a massive 
and very young star. 

6.3. Source n 

As already noted, source n l|Menten fc Reidlll995|) IS an- 
other major power source in the region, in addition to 
source I and BN. Source I and n he on a line running 
NE-SW as shown in fig. ^1 Therefore source n is also lo- 
cated on the axis of the low- v elocity outfl ow in Region B. 
Furthermore IShuping et all t OOi) and iGreenhill etlD 
l,2004b..l ■ see also lBeuther et al. ( 2004) , have made obser- 
vations which suggest the presence of a disk around source 
n with the rotation axis aligned with the low- velocity out- 
flow. Thus source n, which is a possible centre of expan- 
sion of numerous H^O maser s in the vicinity l|Genzel et al.] 
Il98lt iMenten Re iidl ll flfl.'j) . may be the source of or 
may contribute to the low- velocity outflow in Region B. 
However, at present source I woul d seem the clearer candi- 
date following the observations of lGreenhill et all l)2004a|) 
and also the fact that s ource n may have a rel atively low 
luminosity of ~2000 L© llOreenhill et alJ2004b|) compared 
to source I. 

7. Conclusions 

The major conclusions of this study using GriF data for 
the central region of Orion OMCl may be summarised as 
follows. 

(i) 193 distinct bright H2 emission features in the 
v=l-0 S(l) fine have been identifled. 67% of these show 
clear radial velocity structure with well defined features 
of greater than 5 kms~^, indicating that the emission is 
generated by shocks in the medium. 

(ii) There are equal numbers of flows in Peaks 1 
and 2, implying a greater concentration of activity in 
Peak 2 which is less than half the size of Peak 1. Flows 
are found to be grouped much more tightly in Peak 2 
than in Peak 1. Blue shifted flows are somewhat favoured 
over red-shifted, presumably due to dust obscuration. No 
velocities exceed 35-40 kms~^, consistent with the upper 
limit for the propagation of C-type shocks in dense media. 

(in) In a distinct zone. Region B, lying south of 
Peak 1 and west of peak 2, all flows are blue-shifted, with 
an average velocity of — 18±8 kms~^. 

(iv) Weak radial velocity features (WRVs) with ra- 
dial velocities <5 kms^^ are strongly concentrated in 
Peak 1, and a band of such structures has been identified 
in a NE-SW direction across Peak 1. In this band, 2/3 
of the structures show a broken-up morphology, not seen 
elsewhere, which may arise from a major photodissoci- 
ation contribution to the emission through geometrical 
projection effects. This band of structures may also be a 



further manifestation of an outburst from the BN-IRc2 
region about 1000 years ago which created the bullets or 
fingers lying further to the NW. 

(v) Spatial correlation between maximum brightness 
of II2 emission and maximum velocity reveals deviations 
typically of <0.6". These deviations may be used to esti- 
mate the orientation of flows, or shocks, in the plane of 
the sky - but not their absolute direction, that is, whether 
they are travelling north or south, say. Orientations 
are found to be random in Peak 2 and with some clear 
preference for a SE-NW direction in Peak 1, but show a 
strong preference for a NE-SW orientation in Region B. 

(vi) In Paper I it was noted that potential low mass 
star forming zones, identified through the presence of 
energetic flows in OMCl, showed higher densities and 
correspondingly higher energy flux, and are altogether 
less simply structured than more conventional Class 
or 1 low mass isolated star formation regions. Here we 
identify a large number of zones which show similar 
characteristics to those discussed in Paper 1. Data 
point in particular to the close physical association of 
numerous strong H2 emission features. It would seem 
likely that these associations represent sites of multi- 
ple star formation characterised by episodic outbursts 
of activity creating several outflows in different directions. 

(vii) The most striking feature of the results is the 
blue-shifted outflow in Region B. This is the infrared 
counterpart to the outflow associated with the very young 
O-star, source I, a radio source with associated SiO 
and H2O masers and non-maser emission in numerous 
molecules. Radio data show the presence of a blue-shifted 
flow, part of a general expansion with an average ve- 
locity of 18 ±6 kms~^, in agreement with our value of 
-18±8 kms-i. 

(via) There is a remarkable feature in the H2 emis- 
sion data with strong velocity components. On the basis 
of a VLT image, this is suggestive of a region in which 
a flow has burst through the material of OMCl from 
behind, originating in source I and sending shock waves 
through the medium. The axis of the radio outflow from 
source I passes through this object. These and other data 
lead us to conclude that source I is a major progenitor 
of activity in OMCl and in particular of the outflow in 
region B. 

(ix) There remain several extraordinary objects in 
the field whose nature we do not comprehend. For 
example there is a structure, 3" x 3", Table IX!2l object 
2-32 and 2-34 in Peak 2, which is the only clump of 
emission in the entire field with v;.,^ consistently aroun d 
>40 kms"^ (shown in Fig. 4 of ICustafsson et al.ll200(tl . 
The origin of such an energetic object is unknown. 

In conclusion, the GriF data with its associated 
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high spatial and velocity resolution provide information 

which opens a new window on OMCl and presents a 
new perspective on the star-forming activity for both 
massive and low mass stars. Data have very recently been 
obtained for OMCl using the VLT with adaptive optics 
and Fabry-Perot interferometry in the v=l-0 S(l), S(0) 
and the v=2-l S(l) lines. Shock structures seen in these 
3 lines will constrain the physical conditions strongly in 
star forming zones in OMCl. The chief lack at present is 
of shock models which contain both detailed chemistry, 
essential for C-type shocks, and also the geometry, for 
example of bow shocks. Models are under development 
and will be necessary to settle such basic questions 
as the relative physical disposition of the location of 
the maximum brightness and maximum velocity in H2 
emission features. 
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Table A.l. Features identified in OMCl in Peak 1. Data are sorted by the east-coordinate of the maximum in brightness. 
Column 1 labels each feature and columns 2—5 show the position in the plane of the sky of the maxima in brightness and velocity 
for each of the features. Positions are measured in arcseconds east and north of the star TCC0016 (05*' 35™ 14= 91, -05°22'39'.'31 
(J2000) ). Where no coordinates are given for a maximum in velocity, no well-defined peak in the velocity was found nearby. 
Columns 6 shows the velocity of the gas surrounding each feature expressed as visr- Column 7 shows the maximum velocity of 
the flow relative to this local background velocity (see sect. l3.l|l . Column 8 contains the maximum brightness of each feature per 
400 s exposure. Column 9 shows the position angle for each flow, measured east of north (i.e. counter-clockwise), see sect. 15.21 A 
dash indicates that the object is a weak radial velocity region (WRV) and 'no p. a.' indicates a displacement between maximum 
velocity and brightness of <0.25". A star (*) marks flows where the morphology of the emission clearly indicates the absolute 
direction of the flow (see sect. 15. ll : these are only found in Table IX!3l Column 10 gives the displacement in the plane of the 
sky between the observed maxima in brightness and velocity for each feature. A summary of the data can be seen in Fig. |1| 



Flows in Peak 1 





Coordinates of maxima 


velocity of 


flow 


maximum 


position 


displace- 


No. 


Brightness 


Velocity 


surroundings 


velocity 


brightness 


angle 


ment 




E 


N 


E 


N 


kms~^ 


kms^^ 


counts 


degrees 


arcsec 


1-1 


-3.61 


43.58 


— 


— 


17 


- 


204 


- 


- 


1-2 


-7.14 


43.68 


-7.11 


43.12 


16 


-6 


207 


356 


0.6 


1-3 


-7.77 


22.61 


-7.60 


22.54 


8 


-14 


421 


no p. a. 


0.2 


1-4 


-8.51 


46.87 


-8.47 


47.39 


20 


9 


398 


184 


0.5 


1-5 


-9.48 


47.28 


— 


— 


20 


- 


368 


- 


- 


1-6 


-9.52 


27.09 


— 


— 


8 


- 


302 


- 


- 


1-7 


-9.90 


44.94 


— 


— 


21 


- 


302 


- 


- 


1-8 


-10.22 


42.46 


-10.05 


41.90 


21 


7 


264 


343 


0.6 


1-9 


-10.26 


46.38 


- 


— 


20 


- 


366 


- 


- 


1-10 


-12.18 


29.89 


-13.06 


30.73 


13 


-16 


227 


134 


1.2 


1-11 


-12.36 


36.85 


- 


— 


16 


- 


224 


- 


- 


1-12 


-12.39 


28.59 


-11.94 


28.77 


13 


19 


196 


249 


0.5 


1-13 


-13.16 


41.16 


— 


— 


17 


- 


290 


- 


- 


1-14 


-13.61 


29.05 


— 


— 


13 


- 


229 


- 


- 


1-15 


-13.86 


22.61 


-13.48 


22.43 


6 


-28 


187 


294 


0.4 


1-16 


-14.18 


40.95 


— 


— 


18 


- 


328 


- 


- 


1-17 


-14.91 


46.72 


_ 


_ 


22 


- 


259 


- 


- 


1-18 


-15.19 


42.98 


-15.54 


42.84 


14 


-9 


412 


68 


0.4 


1-19 


-15.30 


38.53 






18 


- 


323 


- 


- 


1-20 


-15.75 


32.80 


-16.27 


32.72 


6 


-7 


308 


82 


0.5 


1-21 


-15.86 


30.80 


-15.54 


29.47 


10 


13 


324 


347 


1.4 


1-22 


-16.49 


37.66 






14 




271 






1-23 


-17.50 


43.82 






12 




426 






1-24 


-17.67 


21.25 


-17.53 


21.52 


11 


-10 


236 


207 


0.3 


1-25 


-17.75 


36.51 






11 




410 






1-26 


-18.48 


35.49 


-18.41 


34.76 


9 


-9 


543 


355 


0.7 


1-27 


-18.51 


32.72 






8 




445 






1-28 


-18.59 


23.34 


-17.64 


23.76 


10 


-29 


227 


246 


1.0 


1-29 


-18.69 


42.70 






14 




409 






1-30 


-18.76 


30.38 


-18.83 


29.75 


4 


-8 


362 


6 


0.6 


1-31 


-19.25 


21.21 






12 




304 






1-32 


-19.99 


25.93 






9 




257 






1-33 


-20.37 


28.67 


-20.12 


28.45 


8 


-40 


316 


311 


0.3 


1-34 


-20.58 


35.98 


-20.26 


35.00 


13 


10 


483 


342 


1.0 


1-35 


-20.61 


42.24 






16 




286 






1-36 


-20.68 


30.14 


-20.97 


29.30 


8 


-17 


536 


18 


0.9 


1-37 


-21.00 


32.55 






12 




439 






1-38 


-21.21 


33.99 


-21.14 


33.95 


12 


-8 


439 


no p. a. 


0.1 


1-39 


-21.49 


37.28 


-21.07 


37.45 


14 


-9 


480 


247 


0.5 


1-40 


-21.84 


38.64 


-21.84 


39.87 


15 


-7 


401 


180 


1.2 
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A.l. continued. 





Coordinates of maxima 


velocity of 


flow 


maximum 


position 


displace- 


No. 


Bright 


ness 


Velocity 


surroundings 


velocity 


brightness 


angle 


ment 




E 


N 


E 


N 


kms 


kms 


counts 


degrees 


arcsec 


1-41 


-21.88 


34.51 


— 


— 


12 


— 


416 


— 


— 


1-42 


-22.01 


21.74 






7 


— 


223 


— 


— 


1-43 


-22.23 


24.29 


-21.88 


24.01 


11 


-9 


357 


309 


0.4 


1-44 


-22.43 


30.24 


— 


— 


10 


— 


389 


— 


— 


1-45 


-22.50 


18.48 


— 


— 


7 


— 


471 




— 


1-46 


-22.61 


19.92 


— 


— 


6 




295 






1-47 


-23.80 


35.21 


-23.91 


35.14 


12 


15 


332 


'Hi p. a. 


0.1 


1-48 


-24.11 


37.56 


-23.94 


36.65 


13 


-8 


329 


349 


0.9 


1-49 


-24.33 


19.67 


-23.98 


19.57 


10 


-14 


413 


287 


0.4 


1-50 


-24.57 


21.03 


-24.33 


20.79 


12 


13 


271 


315 


0.3 


1-51 


-24.82 


28.21 


— 


— 


16 


~ 


255 


~ 


~ 


1-52 


-25.48 


29.92 


-25.55 


29.19 


14 


9 


282 


5 


0.7 


1-53 


-25.66 


18.20 


— 


— 


8 


— 


278 


— 


— 


1-54 


-26.04 


35.53 


-25.76 


35.35 


10 


-31 


349 


302 


0.3 


1-55 


-26.53 


24.40 


— 


— 


17 


— 


301 


— 


— 


1-56 


-26.60 


20.79 


-26.28 


20.86 


8 


-5 


270 


257 


0.3 


1-57 


-26.67 


32.13 






14 




378 






1-58 


-27.34 


33.32 


-26.74 


33.39 


12 


-20 


394 


263 


0.6 


1-59 


-27.55 


20.06 


— 


— 


12 


— 


334 


— 


— 


1-60 


-28.17 


34.33 


-28.49 


33.25 


12 


19 


498 


16 


1.1 


1-61 


-29.01 


21.52 


-28.49 


21.49 


8 


-32 


293 


274 


0.5 


1-62 


-29.68 


38.92 


-29.51 


38.81 


19 


-8 


296 


no p. a. 


0.2 


1-63 


-29.68 


23.66 


— 


- 


10 




260 


~ 


~ 


1-64 


-30.73 


25.31 


-31.11 


25.66 


15 


8 


257 


132 


0.5 


1-65 


-30.84 


22.50 


— 


— 


11 


— 


276 


— 


— 


1-66 


-31.50 


20.75 


- 


- 


13 


— 


335 


— 


— 


1-67 


-31.89 


37.62 


— 


— 


26 




235 






1-68 


-31.99 


31.57 


-31.92 


31.11 


17 


-13 


266 


351 


0.5 


1-69 


-32.24 


25.90 


— 


— 


18 


— 


257 


— 


— 


1-70 


-32.69 


24.75 


— 


— 


18 


— 


251 


— 


— 


1-71 


-33.39 


33.74 


— 


— 


28 


— 


241 


— 


— 


1-72 


-33.78 


27.34 


-34.79 


27.23 


20 


15 


268 


84 


1.0 


1-73 


-35.67 


39.24 


-35.35 


38.88 


18 


-11 


270 


318 


0.5 


1-74 


-35.84 


35.07 


— 


— 


27 


— 


438 


— 


— 


1-75 


-36.22 


28.70 


-36.92 


28.63 


23 


-16 


178 


84 


0.7 


1-76 


-36.65 


32.97 


-35.74 


32.65 


22 


-16 


228 


289 


1.0 


1-77 


-37.28 


21.39 






18 


— 


235 


— 


— 


1-78 


-37.38 


24.29 






22 




236 






1-79 


-37.45 


32.97 






24 




263 






1-80 


-37.73 


23.24 






19 




237 






1-81 


-38.88 


32.94 


-39.94 


32.83 


24 


7 


350 


84 


1.1 


1-82 


-39.17 


25.83 






22 




296 






1-83 


-40.60 


22.19 






18 




230 






1-84 


-40.95 


24.29 






16 




178 






1-85 


-44.80 


28.94 






21 




179 






1-86 


-45.29 


35.46 


-45.57 


36.12 


24 


33 


178 


157 


0.7 


1-87 


-49.17 


26.67 






20 




185 
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Table A. 2. Flows in Peak 2. For an explanation of each column, see table IA.l1 



Flows in Peak 2 





Coordinates of maxima 


velocity of 


flow 


maximum 


position 


displace- 


No. 


Brig 


htness 


Velocity 


surroundings 


velocity 


brightness 


angle 


ment 




E 


N 


E 


N 


kms^'^ 


kms~^ 


counts 


degrees 


arcsec 


2-1 


19.39 


-6.51 


19.25 


-6.62 


20 


14 


261 


no p. a. 


0.2 


2-2 


19.00 


9.45 


19.42 


9.41 


34 


9 


233 


275 


0.4 


2-3 


18.59 


-5.43 


18.73 


-5.36 


20 


9 


508 


no p. a. 


0.2 


2-4 


18.38 


3.67 


17.75 


3.71 


34 


-12 


532 


93 


0.6 


2-5 


17.64 


-0.74 


17.40 


-1.23 


24 


-7 


391 


27 


0.5 


2-6 


16.76 


8.72 


16.84 


8.72 


41 


7 


407 


no p. a. 


0.1 


2-7 


16.34 


9.73 


- 


- 


41 


- 


391 


- 


- 


2-8 


16.20 


-6.47 


16.10 


-6.51 


29 


10 


357 


no p. a. 


0.1 


2-9 


16.17 


-0.49 


16.66 


-0.67 


24 


-9 


331 


290 


0.5 


2-10 


15.82 


-1.61 


15.47 


-1.12 


24 


10 


860 


144 


0.6 


2-11 


14.81 


6.58 


- 


- 


34 


- 


287 


- 


- 


2-12 


14.18 


9.80 


13.90 


9.52 


40 


-8 


210 


45 


0.4 


2-13 


13.48 


-1.75 


- 


- 


30 


- 


475 


- 


- 


2-14 


13.37 


-0.46 


13.51 


-0.42 


30 


7 


506 


no p. a. 


0.1 


2-15 


12.91 


8.26 


- 


- 


35 


- 


209 


- 


- 


2-16 


12.36 


0.67 


- 


- 


32 


- 


454 


- 


- 


2-17 


12.22 


-2.66 


12.36 


-3.12 


25 


-6 


244 


343 


0.5 


2-18 


12.15 


2.28 


- 


- 


32 


- 


436 


- 


- 


2-19 


11.65 


3.95 


11.62 


3.88 


32 


-8 


404 


no p. a. 


0.1 


2-20 


10.81 


1.23 


10.57 


1.64 


32 


11 


456 


150 


0.5 


2-21 


10.53 


-0.94 


10.43 


-0.31 


32 


8 


361 


171 


0.6 


2-22 


10.40 


8.75 


10.81 


8.78 


34 


-7 


188 


265 


0.4 


2-23 


10.11 


-6.51 


- 


- 


23 


- 


273 


- 


- 


2-24 


9.66 


1.19 


9.84 


1.01 


31 


7 


349 


no p. a. 


0.2 


2-25 


9.56 


2.17 


9.38 


3.01 


31 


-10 


341 


168 


0.9 


2-26 


8.86 


1.78 


8.78 


1.78 


31 


-12 


370 


no p. a. 


0.1 


2-27 


7.70 


6.30 


8.12 


6.20 


30 


-7 


220 


284 


0.4 


2-28 


7.53 


-3.64 


- 


- 


27 


- 


308 


- 


- 


2-29 


7.07 


-1.05 


7.03 


-1.01 


30 


-6 


335 


no p. a. 


0.0 


2-30 


6.82 


-4.69 


- 


- 


26 


- 


308 


- 


- 


2-31 


6.62 


6.86 


6.51 


6.65 


26 


-11 


182 


no p. a. 


0.2 


2-32 


6.62 


2.62 


6.72 


2.73 


30 


32 


283 


no p. a. 


0.1 


2-33 


6.51 


-3.05 


6.47 


-3.01 


28 


-6 


350 


no p. a. 


0.0 


2-34 


6.30 


1.09 


6.86 


1.33 


30 


26 


310 


246 


0.6 


2-35 


6.30 


-0.56 


6.30 


-0.42 


30 


15 


498 


no p. a. 


0.1 


2-36 


5.95 


7.18 


5.99 


7.24 


26 


-13 


194 


no p. a. 


0.1 


2-37 


5.25 


-3.33 


5.14 


-3.05 


32 


10 


530 


159 


0.3 


2-38 


5.07 


-2.77 


5.07 


-2.24 


30 


9 


323 


180 


0.5 


2-39 


4.80 


-8.23 


4.24 


-7.98 


20 


-7 


205 


114 


0.6 


2-40 


4.62 


10.85 


4.41 


10.95 


20 


-18 


531 


no p. a. 


0.2 


2-41 


3.82 


-2.38 


3.64 


-2.56 


31 


11 


209 


no p. a. 


0.2 


2-42 


3.12 


7.66 


3.01 


7.45 


20 


-8 


305 


no p. a. 


0.2 


2-43 


2.62 


0.94 


2.66 


1.09 


20 


-7 


232 


no p. a. 


0.1 


2-44 


1.92 


1.54 


1.68 


1.58 


20 


-14 


227 


no p. a. 


0.2 


2-45 


1.58 


5.28 


1.78 


4.93 


18 


6 


542 


329 


0.4 


2-46 


1.37 


4.03 


1.92 


3.67 


15 


4 


377 


302 


0.7 


2-47 


0.74 


6.05 


0.77 


6.41 


10 


-13 


372 


186 


0.4 


2-48 


0.14 


4.03 


-0.42 


4.13 


10 


-8 


316 


101 


0.6 


2-49 


-0.07 


8.05 


-0.14 


8.30 


12 


-15 


229 


164 


0.3 


2-50 


-1.05 


9.80 






16 




192 
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Table A. 3. Flows in Region B. For an explanation of each column, see table I A. l1 



Flows in Region B 





Coordinates of maxima 


velocity of 


flow 


maximum 


position 


displace- 


No. 


Bright 


ness 


Velocity 


surroundings 


velocity 


brightness 


angle 


ment 




E 


N 


E 


N 


kms^^ 


kms"'^ 


counts 


degrees 


arcsec 


B-1 


-2.20 


5.60 


-2.13 


5.46 


4 


-14 


240 


no p. a. 


0.2 


B-2 


-3.26 


-12.43 


-3.61 


-12.36 


22 


-9 


174 


101 


0.4 


B-3 


-3.99 


6.12 


-3.61 


6.72 


5 


-22 


348 


213 


0.7 


B-4 


-4.4f 


-6.05 


-4.83 


-5.70 


22 


-8 


192 


130 


0.5 


B-5 


-4.66 


5.70 


-5.57 


5.64 


5 


-34 


432 


86 


0.9 


B-6 


-4.93 


-8.89 


-5.04 


-8.89 


22 


-7 


145 


no p. a. 


0.1 


B-7 


-7.77 


7.80 


-7.95 


7.84 


7 


-28 


208 


no p. a. 


0.2 


B-8 


-8.05 


13.90 


-7.98 


13.69 


14 


-24 


255 


no p. a. 


0.2 


B-9 


-8.05 


7.66 


-7.24 


7.32 


7 


-20 


199 


293 


0.9 


B-10 


-8.19 


7.32 


-7.53 


7.00 


7 


-24 


198 


295 


0.7 


B-11 
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